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1  Introduction 


A  vaxiety  of  problems  in  continuum  mechanics  require  the  solution  X  of  a  linear 
algebraic  equation  of  the  form 

AX +  XA  -  ^(A,  H).  (1.1) 

Here  A,  X,  and  H  are  second-order  tensors  (i.e.,  linear  transformations)  cn  a  two-  or 
three-dimensional  inner  product  space  V,  and  ^(A,H)  is  an  isotropic  function  of  A 
and  H  which  is  linear  in  H. 

For  example,  consider  a  smooth  motion  with  deformation  gradient  F.'  By  the 
polar  decomposition  theorem,  F  has  the  unique  multiphcative  decompositions 

F=RUn.VR,  (1.2) 

where  the  proper  orthogonal  tensor  K.  is  the  rotation  tensor,  and  the  symmetric, 
positive-definite  tensors  U  and  V  are  the  right  and  left  stretch  tensors.  Let  C  and  B 
denote  the  right  and  left  Cauchy-Green  tensors: 

C  =  F^F  =  U2  and  B  -  FF"^  =  .  (1,3) 

Let  the  stretching  tensor  D  and  the  spin  tensor  W  denote  the  symmetric  and  skew 
parts  of  the  velocity  gradient  L: 

L=:D  +  W,  D  =  symL  =  i(L-f  L'T),  W  -  skw  L  =  ^  (L  ~  L^) ,  (1.4) 

•  0 

For  any  tensor  field  A,  let  A  denote  the  material  time  derivative  of  A,  emd  let  A 

denote  the  Jaumann  rate  of  A: 

A~A4AW-WA.  (1.5) 

Then  the  material  time  derivatives  of  the  stretch  tensors  are  related  to  the  material 
time  derivatives  of  the  Cauchy-Gr«*en  tensors  by  the  equations 

UU  +  iJU  -  C  and  VV  4  VV  :=  B  .  ( 1.6) 

The  material  time  derivative  of  the  left  stretch  tensor  is  related  to  the  velocity  gradient 
by  the  equation 

VV  4  VV  VUT  4  ,  (1,7) 

and  the  Jaunianu  rate  of  the  lefi  st fetch  tensor  is  related  to  the  stretching  tensor  by 
the  equation 

VV  4  VV  =  V^D  4  DV'G  (1.8) 

‘We  use  tlie  notation  and  terminology  of  Ttucsdell  and  Noli  [1];  cf  also  Wang  and  llue.sdf;!!  f‘2], 
Gurtin  [3],  and  l>ue,‘jdell  [4] 
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The  material  time  derivative  of  the  right  stretch  tensor  is  related  to  the  tensor 
Du  =  R^DR  by  the  equivalent  equations 

UU  +  UU  =  2UDbU  and  U-^U  +  UU-i  =  2Dr  .  (1.9) 

The  stretching  tensor  is  related  to  the  Jaumauin  rate  of  the  left  Cauchy-Green  tensor 
or  its  inverse  by  the  equivalent  equations 

BD4^DB  =  6  aixd  B-ID  +  DB-’ = -(B->)°.  (1.10) 

The  skew  tensor  fl  =  RR^  is  related  to  the  velocity  gradient  by  the  equation 

vn-f  nv  =  Lv- VL^,  (i.ii) 

and  the  difference  of  W  and  ft  is  related  to  the  stretching  tensor  by  the  equation 

V(W-D)-|-(W-D)V  =  VD-DV.  (1.12) 

The  tensor  equations  (1.6)-(1.12)  have  been  studied  by  various  authors;  cf. 
Leonov  [5],  Sidoroff  [6],  Dienes  [7],  Guo  [8],  Hoger  and  Carlson  [9],  Roger  [10], 
Mehrabadi  and  Nemat-Nasser  [11],  Stickforth  and  Wegener  [12].  and  Guo,  Lehmann 
and  Liang  [13].  These  equations  &ie  of  the  general  form  fl.l)  with  A  =  V,  U,  U“C 
B,  or  B~*,  and  with  $(A,H)  of  the  form 

H,  A’HT-i-  HA^  A^H  +  HA^  AHA,  HA  -  AH^,  AH  -  HA  (1  13) 


In  particular,  for  the  kinematics  applications  discussed  above,  the  coefficient  tensor  A 
in  (1.1)  is  symmetric  and  positive-definite.  These  restrictions  on  A  will  be  assumed 
for  the  present  discussion  only.  They  guarantee  that  a  solution  X  exists  and  is  unique. 
Indeed,  relative  to  any  principal  basis  {e,}  for  A,  the  components  of  X  are  given  by 
the  simple  formula 


X 


•} 


U,  -t  Qj  ' 


(1.14) 


where  a,  is  the  (necessarily  positive)  eigenvalue  of  A  corresponding  to  e,,  and 
are  the  components  cf  ♦(A,H)  relative  to  {e,}.  Observe  that  X  is  symmetric  (resp. 
skew)  iff  ^(A,  H)  is  symmetric  (resp.  skew).  Of  course,  to  actually  compute  X  by 
means  of  (1-14)  we  must  first  determine  the  eigenvalues  and  eigenvectors  of  A. 

For  problems  in  which  the  eigenvalues  and  eigenvectors  of  A  are  not  of  primary 
interest,  it  may  be  more  useful  to  express  X  directly  in  ternas  of  the  tensors  A  and 
H.  Explicit  solutions  of  the  this  type  have  been  derived  by  the  authors  cited  above. 
For  example,  Sidoroff  [6]  and  Guo  [8]  obtained  the  following  solution  of  the  tensor 
equation 

AXfXA  =  H  (1.15) 
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for  the  case  where  H  (and  hence  X)  is  skew  and  dim  V  =  3: 

(/*//*  -  ///a)X  =  {IjI  -  //*)  H  -  (A’H  +  HA^  .  (1.16) 

Here  7^,  77ai  and  IIIx  denote  the  principal  invariants  of  A,  and  the  requirement 
that  A  be  positive-definite  guarantees  that  ~  positive.  Sidoroff  and  Guo 

arrived  at  this  solution  by  first  deriving  a  formula  for  the  ajcial  vector  of  X  in  terms 
of  the  ajciaJ  vector  of  H  and  then  converting  this  intermediate  result  to  its  equivalent 
tensor  form  (1.16).  Hoger  and  Carlson  [9]  obtained  the  following  solution  of  (1.15) 
for  arbitrary  H  when  dim  V  =  3: 

2/77a(/a//a  - //^a)X  =  /aA^HA^ -7a"(A2HA  + AHA2) 

-f  (7a//a  -  III  a)  (  A^H  -f  HA2) 

-f-  +  J77a)  aha  -  7A^7iA(AH  -f  HA) 

+  [/a"7/7a  +  IIxilAllA  -  IIIa)]  H  .  (1.17) 

Since  equations  of  the  form  (1.16)  and  (1.17)  are  said  to  be  displayed  in  direct  no¬ 
tation,  we  will  refer  to  such  equations  as  direct  formulas  for  X  or  direct  solutions 
of  (1.15).  By  a  general  direct  solution  of  (1.15)  we  mean  a  solution,  such  as  (1.17), 
which  is  valid  for  any  tensor  H. 

Although  the  component  formula  (1.14)  is  easily  derived  and  is  independent  of 
the  dimension  of  V,  the  derivation  of  direct  formulas  for  X  is  nontrivial,  and  the 
complexity  of  these  formulas  increases  rapidly  with  the  dimension  of  V.^  For  example, 
when  dim  V  =  2  the  solution  of  (1.15)  for  skew'  H  is^  7aX  =  H,  which  is  substantially 
simpler  than  its  three-dimensional  counterpart  (1.16).  Also,  observe  that  there  is  no 
apparent  simplification  of  the  direct  formula  (1.17)  when  H  is  skew;  in  particular,  it 
is  by  no  means  obvious  that  (1.17)  and  (1.16)  are  equivalent  for  skew  H.  By  utilizing 
Rivlin’s  [17]  identities  for  tensor  polynomials  in  two  variables,  Hoger  and  Carlson  [9] 
were  able  to  convert  (1.17)  to  a  form  which  does  indeed  collapse  to  (1.16)  when  H  is 
skew. 

This  paper  is  devoted  to  the  derivation  and  applications  of  direct  solutions  of 
the  tensor  equation  (1.1)  in  three  dimensions.  Clearly,  for  any  function  $(A,H)  we 
can  obtain  a  direct  solution  of  (1.1)  by  replacing  H  with  #(A,H)  in  (1.17)  or  in 
any  other  general  direct  solution  of  (1.15).  The  resulting  formulas  will  typically  be 
more  complicated  thc:i  the  direct  solution  of  (1.15)  from  which  they  were  obtained, 

^Cf.  the  direct  solutions  of  (1-15)  obtained  by  Smith  [14],  Jameson  [15],  and  Muller  [16].  Their 
formulas  are  valid  for  arbitrary  dimensions,  but  the  complexity  of  these  formulas  is  such  that  they 
would  seem  to  be  useful  only  for  dimV  <  3  or  4.  Also,  compare  Hoger  and  Carlson’s  [9]  solutions  of 
(1.15)  in  two  and  three  dimensions,  and  Mehrabadi  and  Nemat- Nasser’s  [11]  solutions  of  (1.19)  in 
two  and  three  dimensions. 

^This  solution  is  a  special  case  of  the  second  of  two  genera!  direct  solutions  of  (115)  obtained  by 
Hoger  and  Carlson  [9]  in  two  dimensions. 
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although  subsequent  applications  of  the  Cayley-Hamilton  theorem  or  Rivlin’s  [17] 
identities  may  yield  substantial  simplifications  in  some  cases.  One  of  the  goals  of 
this  paper  is  to  develop  methods  which  yield  these  simpler  formulas  more  directly. 
Another  goal  is  to  derive  the  skew  solution  (1.16)  and  other  simple  solutions  of  (1.1) 
for  skew  ^(A,H)  without  resorting  to  intermediate  results  in  terms  of  axial  vectors 
or  to  the  more  complicated  general  direct  solutions. 

The  paper  is  organized  as  follows.  In  Section  2  we  study  the  fourth-order  tensors 
La,  Ma,  and  Na  characterized  by  the  conditions 

X  =  La[H]  AX  +  XA  =  H,  (1.18) 

X  =  Ma[H]  <!=?►  AX-fXA  =  AH-HA,  (1.19) 

X  =  Na[H]  AX -f  XA  =  A^H  -  2AHA -h  HA2 .  (1.20) 

Then  X  is  the  solution  of  the  tensor  equation  (1.1)  iff  X  =  La[^(A,H)].  In  partic¬ 
ular,  Ma[H]  =  La[AH  -  HA]  and  Na[H]  =  LaIA^H  -  2AHA  +  HA2].  Conversely, 

when  $(A,H)  has  one  of  the  forms  in  (1.13),  we  show  that  there  axe  simple  relations 

for  La[^(A,H)]  in  terms  of  Ma[H],  MAfsymH],  or  Na[H].  The  utility  of  these 
relations  is  due  to  the  fact  that  direct  formulae  for  Ma[H]  and  Na[H]  are  simpler 
and  easier  to  derive  than  general  direct  formulas  for  X  =  La[H]  such  as  (1.17).  The 
results  in  Section  2  are  independent  of  the  dimension  of  the  inner  product  space  V. 
Furthermore,  unlike  the  component  formula  (1.14),  these  results  are  valid  for  any 
tensor  A  with  the  property  that  (1.15)  has  a  unique  solution  X  for  any  given  H. 
Such  a  tensor  A  is  necessarily  nonsingulax  but  need  not  be  symmetric  or  definite, 

In  Sections  3-7  we  assume  that  dimV  =  3.  Section  3  contains  various  tensor 
identities  which  will  be  utilized  in  the  sequel.  These  include  Rivlin’s  [17]  identities 
for  tensor  polynomi<ils  in  two  variables  as  well  as  some  new  identities  which  follow 
from  Fiiviin’s.  In  Section  4  we  consider  (1.15)  with  X  and  H  restricted  to  the  set 
T(A)  of  2dl  tensors  K  such  that  tr  (A"K)  =  0  (n  =  0, 1,2),  We  obtain  necessary  and 
sufficient  conditions  for  the  existence  of  a  unique  solution  in  T(A)  (the  possibility  of 
other  solutions  outside  T(A)  is  not  excluded  here),  and  we  derive  direct  formulcis  for 
this  solution.  When  A  is  symmetric  these  formulas  are  valid  for  any  skew  tensor  H; 
in  particulcir,  recover  the  formula  (1.16)  of  Sidoroff  and  Guo.  These  results  do  not 
require  that  A  be  nonsingular;  some  applications  for  which  A  might  be  singular  are 
discussed  below.  Section  4  concludes  with  the  derivation  of  necessary  and  sufficient 
conditions  for  the  existence  of  a  unique  solution  X  of  (1.15)  with  H  unrestricted. 
The  proof  utilizes  the  results  for  the  special  case  where  X  and  H  belong  to  T(A).  In 
Section  5  we  use  the  results  in  Sections  3  and  4  to  derive  direct  formulas  for  Ma[H] 
and  Na[H]  for  arbitrary  H;  these  formulas,  together  with  the  relations  for  La  in 
terms  of  Ma  or  Na  obtained  in  Section  2,  are  in  turn  used  to  derive  direct  formulas 
for  La[H]  which  2ure  vcilid  for  arbitrary  H.  Then  direct  formulas  for  La[^(  A,  H)]  with 
as  in  (1.13)  follow  from  these  results  and  the  identities  in  Section  2.  In  Section  6 
we  derive  equations  (1.6)-(1.12)  and  apply  our  results  to  the  solution  of  these  and 
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related  equations  axising  in  the  kinerr  atics  of  continua.  In  Section  7  we  discuss  some 
additional  kinematic  formulas  which  can  be  obtained  from  various  transformations 
of  the  results  in  Section  6.  Although  some  of  the  algebrciic  and  kinematic  formulas 
derived  in  this  paper  have  been  obtained  previously  by  other  authors,  the  derivations 
given  here  are  new,  and  we  derive  many  new  formul  s  as  well. 

Additional  Applications 

The  general  results  in  Sections  2-5  should  prove  useful  for  other  problems  in 
mechanics  which  lead  to  tensor  equations  of  the  form  (1.1).  Some  of  these  problems 
are  listed  below. 

1.  Direct  formulas  for  the  derivatives  of  the  stretch  and  rotation  tensors  with 
respect  to  the  deformation  gradient:  Here  A  =  U  or  V,  and  $(A,H)  has  some  of 
the  forms  listed  in  (1.13)  as  well  as 

AH  +  H^A,  HA  +  ART,  AH-RTA.  (1.21) 

This  problem  is  the  subject  of  a  follow-up  paper  [18];  cf.  also  Wheeler  [19]  and  Chen 
and  Wheeler  [20]  for  a  different  approach  to  this  problem.  For  a  hyperelastic  material, 
the  results  in  [18]-[20]  also  yield  direct  formulas  for  the  first  Piola-Kirchhoff  stress 
tensor  in  terms  of  the  derivative  of  the  strain  energy  function  with  respect  to  the 
right  stretch  tensor  U;  cf.  also  Roger  [21],  where  this  problem  hcis  been  solved  using 
Roger  and  Carlson’s  [9]  formula  (1.17). 

2.  Direct  formulas  for  one  work-conjugate  stress  tensor  in  terms  of  another,  or  for 
a  work-conjugate  stress  tensor  in  terms  of  the  Cauchy  or  first  Piola-Kirchhoff  stress 
tensors  (cf.  Guo  and  Man  [22]):  Here  A  =  U  or  V,  and  ^(A,H)  has  some  of  the 
forms  listed  in  (1.13)  as  well  as 

m 

AH,  HA,  ^A"’-^HA’-’.  (1.22) 

r=l 

3.  The  kinematics  and  dynamics  of  rigid  bodies  (cf.  Truesdell  [4,  §1.10,  1.13] 
and  Scheidler  [23])  and  p.seudo- rigid  bodies  (cf.  Cohen  and  Muncaster  [24]):  Here 
the  symmetric  tensor  A  is  either  the  current  or  the  referential  Euler  tensor.  If  the 
mass  is  not  confined  to  a  single  plane  then  A  is  positive-definite  (Segner’s  Theorem). 
However,  the  results  in  Theorem  4.1  also  apply  when  the  mass  is  confined  to  a  single 
plane,  in  which  case  A  is  singular. 

4.  Traction  boundary  value  problems  in  finite  elasticity:  For  A  symmetric  and  H 
and  X  skew,  the  tensor  equation  (1.15)  arises  in  connection  with  Signorini’s  expansion 
and  Stopelli’s  theorems;  cf.  Wang  and  Truesdell  [2,  §7.2,  7.4].  For  these  applications 
the  astatic  load  tensor  A  may  be  singular;  the  results  in  Section  4  are  applicable  in 
this  case  provided  that  the  load  system  does  not  possess  an  axis  of  equilibrium. 

5.  Stability  analysis  of  systems  of  ordinary  differential  equations:  The  tensor 
equation  A^Y  -f  YA  =  G  arises  in  the  construction  of  quadratic  Liapunov  functions; 
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cf.  Hahn  [25,  Ch.  4]  and  Gantmacher  [26,  §5.5].  At  the  end  of  Section  2  we  show  how 
general  direct  solutions  of  this  equation  can  be  obtained  from  general  direct  solutions 
of  AX  +  XA  =  H. 


2  The  fourth-order  tensors  L^,  M^,  and 

Let  Lin  denote  the  set  of  all  linear  transformations  (or  second-order  tensors)  on 
the  finite- dimensional  real  inner  product  space  V.  Sym  and  Skw  denote  the  subspaces 
of  Lin  consisting  of  all  symmetric  and  skew  tensors,  respectively.  Psym  denotes  the 
set  of  all  symmetric  and  positive-definite  tensors.  The  identity  tensor  is  denoted  by 
I,  and  A°  :=  I  for  any  tensor  A.  Unless  specified  otherwise.  A,  G,  H,  X,  and  Y 
denote  curbitraxy  tensors.  We  assume  that  Lin  is  equipped  with  the  inner  product 
“  •  ”  defined  in  terms  of  the  trace  function  by  H  •  G  ;=  tr  (H^G),  where  denotes 
the  transpose  of  H. 

By  a  fourth-order  tensor  we  mean  a  linear  transformation  from  Lin  into  Lin.  The 
image  of  H  €  Lin  under  a  fourth-order  tensor  K  is  denoted  by  K[H].  In  this  section 
we  study  the  properties  of  the  fourth-order  tensors  L^.,  Ma)  ^nd  discussed  in  the 
Introduction.  To  facilitate  the  statement  and  proof  of  some  of  these  properties,  we 
introduce  the  fourth-order  tensors  and  Ca: 

Ba[X]:=AX  +  XA  and  Ca[X]  :=  AX  -  XA  .  (2.1) 

It  is  easily  verified  that  Bx  and  Cx  commute;  indeed, 

Cx^A  =  Ba^a  —  Cx2  ■  (2.2) 

Clearly,  Cx  is  singular  for  every  tensor  A.  Let  Lin*  denote  the  set  of  all  A  G  Lin 
such  that  B^  is  nonsingular,  or,  equivalently,  such  that  the  equation  AX  -I-  XA  =  H 
heis  a  unique  solution  X  for  any  given  H.  From  the  discussion  in  the  Introduction  we 
know  that  Psym  C  Lin*.  Necessziry  and  sufficient  conditions  for  A  G  Lin*  in  terms  of 
the  principal  invariants  of  A  or  in  terms  of  the  chtiracteristic  roots  of  A  are  discussed 
in  Section  4.  For  this  section  we  need  only  the  following  elementary  results. 

Proposition  2,1  The  following  conditions  are  equivalent: 

(1)  A  G  Lin*; 

(2)  AT  G  Lin*; 

(3)  QAQ~^  G  Lin*  for  every  nonsingular  tensor  Q; 

(4)  A  is  nonsingular  and  A“*  G  Lin*. 

Proof:  The  equivcilence  of  (1)  and  (2)  follows  from  the  equivalence  of  the  equations 
AX  -1-  XA  =  H  and  A^X^  -f  X^A^  =  HT.  Sim  laxly,  the  equivalence  of  (1)  and  (3) 
follows  from  the  equivalence  of  the  equations  AX  -t-  XA  —  H  and 

(QAQ-^)(QXQ-M  +  (QXQ-‘)(QAQ^)  =  QHQh 
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Now  suppose  that  A,  «ind  hence  A'^,  is  singular.  Then  there  are  nonzero  vectors  u 
and  V  such  that  Au  =  0  and  A’^v  =  0.  Therefore, 

(8*  v]  =  (Au)  0  V  +  u  (gi  (A^'v)  —  0 . 

Since  u  (g)  V  7^  0,  is  singular  and  hence  A  ^  Lin*.  Thus  A  e  Lin*  A  is 
nonsingulaj.  That  A  €  Lin*  =»  A“*  G  Lin*  follows  from  the  equivalence  of  the 
equations  AX  f  XA  =  H  and  A~*X  +  XA“'  =  A~^HA~L  Hence,  (1)  implies  (4). 
Conversely,  if  A~'  G  Lin*  then  A  =  (A“*)~*  G  Lin*.  □ 

Proposition  2.1  shows  that  nonsingulaxity  of  A  is  necessary  for  A  G  Lin*;  as  we 
will  see  ii?  Section  4,  it  is  not  sufficient.  For  the.  remainder  of  this  section  we  assume 
that  A  G  Lin*.  Unless  specified  otherwise,  no  additional  restrictions  are  imposed 
on  A.  We  denote  the  inverse  of  Bx  by  L^: 

La--=(Ba)“'.  (2.3) 

Then  Ba[X]  =  H  iff  X  =  La[H],  which  is  equivalent  to  the  statement  (1.18).  If  I 
denotes  the  fourth-order  identity  tensor,  then 

LaBa  =  BaLa  =  I ,  (2.4) 

which  is  equivalent  to  the  relations 

LA(AH-f  HA]  =  ALa[H]  +  La[H]A  =  H.  (2.5) 

If  A  is  symmetric,  then  from  (2.1)i  it  folic-- s  that  Ba,  and  hence  its  inverse  La,  maps 
symmetric  tensors  to  symmetric  tensors  ajid  skew  tensors  to  skew  tensors. 

The  fourth-order  tensor  Ma  is  defined  by 

Ma  :=  LaCa  =  CaLa  ,  (2.6) 

where  (2.6)2  follows  by  multiplying  (2.2)j  on  the  left  and  right  by  La  and  then  using 
(2.4).  From  (2.1)2  if  follows  that  (2.6)  is  equivalent  to  the  relations 

Ma(H|  =  LaIAH  -  HA]  =  ALaIH]  -  La|H]A  .  (2.7) 

By  replacing  H  with  AH  -  HA  in  (1.18),  we  see  that  (2.7)i  is  equivalent  to  the 
statement  (1.19).  Prom  (2.6),  (2.2),  and  (2.4),  we  obtain  the  relations 

MaBa  =  BaMa  =  Ca  =  LaCa3  =  Ca?La  ■  (2.8) 

By  (2.1)  we  see  that  (2.8)  is  equivalent  to  the  relations 

Ma[AH-^HA]  =  AMa[H]  T  Ma[H]A  =  AH  -  HA 

=  LaIA^H  -  HA^l  =  AH.a[H]  -  La[H]A^  .  (2.9) 
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If  A  is  symmetric,  then  AH  -  HA  is  symmetric  (resp.  skew)  if  H  is  skew  (resp. 
symmetric).  Hence,  by  (2.7)i  and  the  coimnent  following  (2.5),  we  see  that  if  A  is 
symmetric  then  maps  symmetric  tensors  to  skew  tensors  and  skew  tensors  to 
symmetric  tensors. 

The  fourth-order  tensor  is  defined  by 

Na  :=  MaCa  =  CaMa  =  La  (C^)'  -  (Ca)'^  ,  (2.10) 

where  (2.10)2-4  follow  from  (2.6).  Since 

(Ca)' [H]  =  A^H  -  2AHA  +  HA2  ,  (2.11) 

(2.10)  is  equivalent  to  the  relations 

Na[H]  =  Ma[AH  -  HA]  -  AMa[H]  -  Ma[H]A 

=  LA[A2H-2AHA-f  HA^j 

=  A^LaIHI-T  La[H]A2-2ALa[H]A.  (2.12) 

By  replacing  H  with  A^H  -  2AHA  -|-  HA^  in  (l.lSj,  we  see  that  (2.12)3  is  equivalent 
to  the  statement  (1.20). 

Proposition  2.2  Let  denote  the  set  consisting  of  the  fourth-order  tensors  intro¬ 
duced  above: 

5a:={Ba.Ca.L*,M*,N*).  (2.13) 

Then  any  two  tensors  in  5a  commute,  and  each  tensor  Ka  €  5a  has  the  following 
properties: 

Ka[-A-iHA2]  =  AiKa[H]A2,  */ Ai  and  A2  commute  with  A,  (2.14) 
Kqaq-‘[QHQ-‘]  =  QKa[H]Q for  any  nonsingular  tensor  Q.  (2.15) 

In  particular,  Ka[H]  is  an  isotropic  function  of  A  and  H  which  is  linear  in  H,  and 

Ka[A"*HA"]  =  A'"Ka[H]A",  for  any  integers  m  and  n.  (2.16) 

Proof:  That  5a  is  commutative  follows  from  (2.2),  (2.4),  (2.6),  (2.8),  and  (2.10). 
The  easiest  way  to  establish  the  other  properties  is  to  first  prove  them  for  the  ten¬ 
sors  Ba  and  Ca,  and  then  use  the  definitions  (2.3),  (2.6),  and  (2.10)  to  prove  the 
corresponding  results  for  La,  Ma,  and  ^a-  We  prove  (2.14)  and  leave  the  proof 
of  (2.15)  to  the  reader.  It  is  easily  seen  that  (2.14)  holds  for  Ka  =  Ba  or  Ca- 
To  prove  (2.14)  for  Ka  =  I-a,  apply  La  to  BaIAiXAj]  =  AiBa[X]A2  to  obtain 
AjXAj  =  La[AiBa[X]A2],  and  then  set  X  =  La[H].  To  prove  (2.14)  for  Ka  =  Ma, 
use  the  fact  that  (2.14)  holds  for  Ka  —  Cj^  and  La^ 

Ma[A,HA2]  =  La[Ca[A2HA2]]  =  La[A,Ca[H]A2] 

=  .A,La[Ca[H]]A2  =  A,Ma[H]A2. 
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The  proof  of  (2  14)  for  is  similar.  □ 

The  equations  (2.7)  give  two  expressions  for  Myv[H]  in  terms  of  L^,  A,  and  H. 
Alternate  expressions  are 

Ma[H]  =  H-2La[HA]  =  H~2La[H]A 

=  2La[AH)-H  =  2ALa[H]-H.  (2.17) 


The  relations  (2.17)i  3  follow  from  (2.7)i  and  (2.5);  for  exajnplo, 

Ma[H]  =  La[AH-HA]  =  Lx[AH4  HA  -2HA] 

=  L^f AH  4  HA]  -  2  La[HA]  =  H  -  2 Lx[HA] . 

Then  (2.17)2,4  follow  from  (2.17)^  3  and  (2.16)  with 

Proposition  2.3  The  fourth- order  tensor  can  be  expressed  in  terms  of  the  fourth- 
order  tensor  by  the  formulas 

U|H1  =  J(H  -  Mx[Hl)A-'  =  1A-‘(H  +  .  (2.18) 

Similarly,  can  be  expressed  in  terms  o/N^  by  the  formulas 

U(H1  =  l(A-'H  +  HA-'-A-'N*|H)A-') 

=  1A-‘(AH  +  HA-Nji[H1)A-'.  (2.19) 

Proof:  (2.18)  follows  from  (2.17).  Then  from  (2.18)  we  have 

L*[H|  =  i(A-‘H  +  HA->  +  A-'M*|H]  -  [HJA  -), 

and  from  (2.12)  we  have 

A-iMa[H] -Ma[H]A-‘  =  A'1(Ma[H]A- AMx[H])A-> 

=  A-M-Na[H])A-'.  n 

Now  suppose  that  by  some  means  we  have  obtained  a  direct  solution  of  the  tensor 
equation  (1.19),  or,  equivalently,  a  direct  formula  for  Mx[H],  which  is  valid  for  any 
tensor  H.  As  we  will  see  in  Section  5,  such  formulas  are  relatively  simple  and  easily 
derived  when  dimV  =  3.  Then  direct  formulas  for  [H]  which  are  valid  for  any 
tensor  H  follow  from  (2.18).'*  Alternatively,  we  can  use  the  relations  (2.12)i  2  a.nd 
the  direct  forrr.ulas  for  M^^fH]  to  obteiin  direct  formulas  for  N^^fH],  and  then  use 
(2.19)  to  obtain  direct  formulas  for  L^[H]  v.'hich  are  veJid  for  any  H.  In  any  of  these 
formulas  we  cam,  of  course,  replace  the  A"^  terms  by  a  polynomiad  in  A  via  the 

*For  the  tensor  equation  A^X  +  XA  =  H  (H  syitnmetric),  Barn  ttt  and  Storey  [27]  obtained  a 
relation  analogous  to  (2.18)i;  cf.  their  equations  (1-2),  (2.1),  and  (2.2).  Their  relation  is  equivalent 
to  (2.18)i  when  both  H  and  A  are  symmetric.  They  did  not  obtain  direct  solutions. 
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Cayley-Hainiltou  thf*orein.  Theso  techniques  'A-ill  he  u.se<l  in  S>'(  tioii  to 
direct  formulas  for  Lj^[H]  in  three  dimensions. 

By  (1.18)  with  H  — »  ♦(A,H),  the  unique  solution  X  of  the  tensor  equation 
AX  +  XA  =  ♦(A,H)  is  given  by  X  =  Lx(^{A.H)].  If  ♦(A.H)  is  an  isotropic 
function  of  A  and  H  which  is  linear  in  H  (in  particular,  if  ♦(A,H)  has  one  of  the 
forms  in  (1.13),  (1-21),  or  (1.22)),  then  by  (2.15)  it  follows  that  X  is  al.so  an  isotropic 
function  of  A  and  H  which  is  linear  in  H.  When  ♦(  A,  H)  -  AH  —  HA,  this  solution 
can  also  be  written  as  X  =  Mj^[H].  As  we  will  see  in  Section  5,  for  arbitrary  H  the 
direct  formulas  for  Mx[H]  are  much  simpler  than  the  direct  formulas  for  i_j^[Hj,  which 
should  not  be  too  surprising  in  view  of  the  method  described  above  for  generating 
the  latter  formulas.  For  other  functions  ♦(A,H)  in  the  list  (1.13),  the  existence 
of  relatively  simple  direct  formulas  for  X  is  due  to  the  fact  that  there  are  simple 
expressions  for  Lx[^(A,H)]  in  terms  of  Mx{H],  MxfsymH],  or  Nx[H].  We  derive 
these  identities  below.  Similar  results  hold  for  ^(A,  H)  of  the  form  (1.21)  and  (1.22); 
cf.  Scheidler  [18]. 

Consider  the  case  $(A,H)  AHA,  i.e.,  the  tensor  equation 

AX  +  XA  =  AHA.  (2.20) 

Alternate  expressions  for  the  solution  X  =  L^fAHA]  axe  given  by  the  following 
identities: 


U|AHA|  =  AL*[H1A  =  L*-,|H1 

=  1A(H  -  MxIHI)  =  1(H  +  M*[H))A 
=  J(AH  +  HA-N>,[HJ).  (2.21) 

(2.21) i  follows  from  (2.16);  then  (2.21)3  .4  5  follow  from  (2.18)  and  (2.19)2-  To  obtain 

(2.21) 2,  observe  ♦^hat  (2.20)  is  equivalent  to  the  tensor  equation 

A-‘X  +  XA-'=H,  (2.22) 

and  that  the  solution  of  (2.22)  is  X  =  Lx-i[H].  For  A  G  Psym  and  K  G  Sym.,  rela¬ 
tions  equivalent  to  some  of  those  in  (2.21)  were  observed  by  Mehrabadi  and  Nemat- 
Ncisser  [11]®  and  Cohen  and  Muncaster  [24,  Ch.  6].®  Our  derivations  above  and  in 

®Iii  their  analysis  of  the  tensor  equation  (1.9)2  for  U,  Mehrabadi  and  Nemat-Nasser  obtained 
relations  which,  for  symmetric  A  and  H,  are  equivalent  to  the  relations  Lx->[f^]  =  ^(H-f  Ma[H])A 
and  Lx->[H]  =  ^(AH-t-HA  — Na[H]);  cf.  equations  (8.8),  (8.12),  (8.13),  and  (8.16)  in  their  paper. 
They  also  derived  a  direct  formula  for  Ma[H]  in  three  dimensions  and  used  this  formula,  together 
with  the  latter  of  the  two  relations  above,  to  obtain  a  direct  solution  of  (1.9)2;  cf.  (6.12)i  in  this 
paper,  which  we  will  obtain  by  essentially  the  same  technique.  However,  our  derivation  of  direct 
formulas  for  Ma[H]  in  Section  5  differs  substantially  from  the  method  used  in  [11]. 

®Cohen  and  Munc.Tster  considered  a  tensor  equation  of  the  form  A“*X-f-XA"‘  -f-c(trX)A“^  = 
G,  where  A  G  Psym  is  the  referential  Euler  tensor  and  G  is  synunetric.  This  equation  arises  in 
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the  Df^xt  paragraph  differ  from  theirs;  iit  particular,  we  do  not  rely  on  the  symmetry 
of  A  or  H. 

Compared  with  (2.21)2,  the  formula  for  in  terms  of  is  much  simpler: 

Ma-.  --Ma.  (2.23) 

Indeed,  from  (1.19)  we  see  that  X  =  Ma[H]  iff  A~’X+XA~*  =  A"*(  — H)  — (— H)A~* 
iff  X  =  Ma->(-H].  Note  that  some  of  the  identities  in  (2.21)  can  also  be  obtained 
by  replacing  A  with  A“*  in  (2.18)  and  (2.19)  amd  then  using  (2.23)  and  (2.12). 

For  the  case  4p(A,H)  =  A^H  +  HA^,  we  have  the  identities 

La(A^H  +  HA2]  =  A2La[H]  +  La[H]A2 

=  2  La[H]A2  +  ah  -  ha  =  2A2La[H]  -  AH  +  HA 
=  AH  -  Ma[H]A  =  HA  +  AMa[H] 

=  |(AH-1-HA  +  Na[H]).  (2.24) 

(2.24) i  follows  from  (2.16);  then  (2.24)2,3  follow  from  (2.9).  (2.24)4,5  follow  from 

(2.24) 2,3  and  (2.13).  Finally,  (2.24)6  follows  from  (2.12)4  (2-21)5,  or  from  (2.24)4,5 

and  (2.12).  Next,  consider  the  cause  $(A,  H)  =  A^H'^+HA^.  By  applying  La  to  the 
identity 

A^HT+  HA2  =  A2(symH)  +  (symH)A2  -f-  A2(-skwH)  -  (-skwH)A2 

and  then  using  (2.9)  with  H  — skw  H,  we  obtaiin  the  identity 

LaIA^HT+HA^)  =  LA[A2(symH)  +  (symH)A2] 

+  (skwH)A  —  A(skwH) .  (2.25) 

Alternate  expressions  for  LAfA^H'^d-  HA^]  follow  from  this  aind  (2.24)  with  H  — » 
symH.  In  particular,  we  have 

La(A2HT+  HA*]  +  A(skw  H)  -  (skw  H)A 
=  A(symH)  —  MAfsymHjA  =  (symH)A  +  AMA[symH] 

=  |(A(symH)  +  (symH)A  +  NA[symH]) .  (2.26) 

Finally,  for  the  case  #(A,H)  =  HA  —  AH'^,  we  have  the  identity 

La[HA  -  AHT]  =  skw  H  -  Ma [sym H] .  (2.27) 

the  analysis  of  gyroscopic  motions  of  pseudo-rigid  elastic  bodies.  On  multiplying  this  equation  by 
A  and  taking  the  trace  of  the  result  we  may  solve  for  tr  X  and  reduce  the  original  equation  to  the 
form  (2.22)  with  H  =  G  —  c(2 -|- 3c)~Hr  (GA).  Then  their  equations  (6.3.12),  (6.3.14),  and  (6.3.16) 
are  equivalent  to  our  relation  La-*[H]  =  ^(^  +  ^^a[H])A.  They  did  not  obtain  direct  formulas  for 
Ma[H]  or  La-.[H]. 
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This  follows  from  the  identity 

HA  —  AH^  =  —  A(sym H)  +  (symH)A  -f  A(skw  H)  +  (skw  H)A  , 

(2.5)  with  H  —*  skwH,  and  (2.7)i  with  H  — »  symH. 

We  conclude  this  section  with  a  result  mentioned  in  the  introduc  cion  in  connection 
with  Liapunov  functions  for  systems  of  dilferential  equations. 

Proposition  2.4  The  tensor  equation  AJY  -f-  YA  =  G  has  a  unique  solution  Y  for 
any  given  G  iff  A  E  Lin*.  If  A  E  Lin*,  then  X  =  X)m,n  is  a  general 

direct  solution  o/AX  +  XA  -  H  iff  Y  =  !Cm,n  Cm,„(A'^)”^GA"  is  a  general  direct 
solution  of  AJY  +  YA  =  G.  Here  m  and  n  are  integers,  the  sums  are  assumed  to  be 
finite,  and  the  coefficients  c^  „  may  depend  on  A  but  not  on  H  or  G. 

Proof:  We  use  the  well-known  fact  that  A^  and  A  are  similar,  i.e.,  there  is  a 
nonsingular  tensor  S  such  that  AI  =  SAS~L  Let  Y  =  SX  and  G  =  SH.  Then  the 
equations  AX  -I-  XA  =  H  and  A^  Y  -f  YA  =  G  are  equivalent,  and  the  results  of  the 
proposition  follow.  □ 

Since  AX  -H  XA  =  H  iff  X  =  La[H],  Proposition  2.4  can  be  used  to  transform 
the  general  direct  formulas  (5.6)  and  (5.17)-(5.20)  for  L^[H]  in  three  dimensions  into 
general  direct  solutions  Y  of  A'^Y  -f  YA  =  G.  This  proposition  cannot  be  applied 
to  the  general  direct  formulas  (5.22),  (5.28),  and  (5.33)  since  some  of  the  coefficients 
in  these  formulas  depend  on 


3  Some  tensor  identities  in  three  dimensions 

The  derivations  of  the  results  in  Sections  4-7  utilize  various  identities  involving 
one  or  two  second-order  tensors  and  their  principal  invariants.  For  convenience  we 
have  collected  most  of  these  identities  in  the  present  section.  With  the  exception  of 
some  comments  at  the  end  of  Section  4,  for  the  remainder  of  this  paper  we  assume  that 
the  underlying  inner  product  space  V  is  three-dimensional.  Unless  specified  otherwise, 
the  tensors  A  and  H  are  arbitrary. 

The  principal  invariants  of  A  are  denoted  by  /^,  Z/^,  emd  IIIx,  and  its  charac¬ 
teristic  roots  (in  the  complex  field)  are  denoted  by  01,02,03.  Then 

3 

det  (il  -  A)  =  x^~  Ij,x^  -I-  lUx  -  =  n(^  -  “.)  >  (31) 

i=l 

^If  in  Proposition  2.4  we  allowed  the  coefficients  Cm,n  in  the  formula  for  X  to  depend  on  H, 
say  Cm,n  =  Cm,n(A,H),  then  since  H  =  S~‘G  it  follows  that  the  coefficients  in  the  corresponding 
formula  for  Y  would  depend  not  only  on  A  and  G  but  also  on  the  tensor  S  in  the  similarity 
transformation  A^  =  SAS~*. 
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where  (3.1)i  holds  for  any  real  number  i,  (3.1)2  bolds  for  any  complex  x,  and 


/x  =  tr  A  =  a,  +  02  +  03  ,  (3.2) 

lU  =  2  -  ^A?)  =  ai«2  +  a2«3  +  aaOi  ,  (3.3) 

///x  =  det  A  =0,0203.  (3.4) 

The  characteristic  root  o,  is  real  iff  o^  is  an  eigenvalue  of  A.  The  second  and  third 
moments^  of  A  are 


U,  =  I^-2IU  =  a-\^a\  +  a.l,  (3.5) 

/a-  =  /a"  -  3/aWa  +  3///a  =  oj  +  aj  +  (3.6) 

The  Cayley- Hamilton  theorem  implies  that 

A=  =  /*A»-H*A  +  W/aI.  (3.7) 

The  adjoint  of  A  is  the  tensor 

adj  A  :=  A^  -  UA  +  II jJ  =  UIj,A~^  ,  (3.8) 


where  the  expression  on  the  right  is  valid  only  when  A  is  nonsingular. 

Since  the  expression  7x1- A  occurs  frequently  in  the  sequel,  we  introduce  a  special 
symbol  for  it: 

A:=/xI-A.  (3.9) 

Then 

7^  =  2/x  and  A  =  |7;^I  —  A  .  (3.10) 

Expressions  for  the  determinant  of  A  in  terms  of  A  are 

Ilh  =  ^aWa  -  Wa  =  -  ^A-) 

=  (o,  +  02)(a2  -f  03)(03  -f-  Oi).  (3.11) 

(3.11)i  3  follow  from  (3.1)  with  x  =  7x,  and  (3.11)2  follows  from  (3.6),.  (3.11)3  ^Iso 
follows  from  the  fact  that  the  characteristic  roots  of  A  are®  a,  -f  03,  02  +  03,  03-1-01. 
If  A  is  nonsingular,  then  from  (3.8)2,  (3.3),,  and  (3.11)i,  we  obtain*® 

77x  =  777x7x-.  and  777^  =  777x(7x7x-i  -  1) .  (3.12) 

®Cf.  Ericksen  (28,  §38]. 

®For  any  scalar  polynomial  p(i),  the  characteristic  roots  of  the  corresponding  tensor  polynomial 
p(A)  are  p(ai);  cf.  Theorem  3,  p.  84,  in  [29].  By  choosing  p(x)  =  /x  —  *  it  follows  that  the 
characteristic  roots  of  A  are  Ix  —  =  Oj  -1-  ot  for  distinct  i,J,  k. 

*°Cf.  Stickforth  [30]. 
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Expressions  for  the  second  principal  invariant  c*'  A  in  terms  of  A  arc 

//*  =  /*’  +  //*=  i  (3V -/*.)=/*.  + 3//* 

=  (Oi  +  a2)(a2  4-03)  +  (oj -♦-a3)(a3  +  a,)  +  (03  +  aj)(a,  +  “2)  •  (313) 

By  replacing  A  with  A  in  (3.8)  and  using  (3.9),  (3.10)i,  and  (3.13)i,  we  obtain  the 
following  expressions  for  the  adjoint  of  A: 

adj  A  =  A*  +  IIjJ  =  ///^A-*  ,  (3.14) 

where  the  expression  on  the  right  is  vaJid  only  when  A  is  nonsinguiax.”  Let 

A*  :=  I^-IU  =  lA>+IlA  =  i{lA+lA?)  =  iIi,- 

=  i((a, +a2)^  +  (a2  +  a3)^  +  (a3  +  ai)^]-  (3.15) 

This  invariant  appears  in  the  direct  solution  (1.16)  of  Sidoroff  and  Guo  and  also  in 
several  other  formulas  in  the  sequel.  Observe  that  2A_4  is  the  second  moment  of  A. 
The  following  identities  are  due  to  Rivlin  [17]: 

A2HA2  =  //xAHA  - ///x(AH  +  HA) 

+  .^A*hA^  +  Oa^hA  +  >  (3.16) 

A^HA  +  AHA^  =  /aAHA-/7/aH 

+  /*hA» A +  (3.17) 

A^H  +  HA^  +  AHA  =  /*(AH  +  HA)  -  II aH 

+  /„A7  +  Qi,'„A  +  <l„I,  (3.18) 

where 


‘^A.H 

(3.19) 

“a,h 

=  ^A»H  ~  ^A^AH  =  ~-^AaH  7 

(3.20) 

Q'aIh 

=  ^AH  ~  ^A^H  =  ~^AH  7 

(3.21) 

=  ^A>H  ~  ^A^AH  +  IIa^H  =  ^(adjA)H  • 

(3.22) 

The  first  expressions  in  (3.19)-(3.22)  are  the  ones  given  by  Rivlin  [17];  the  second 
expressions  follow  easily  from  these  and  (3.7)-(3.9).  From  (3.9),  (3.15)i,  and  (3.18), 
we  obtain  the  identities 

AHA  =  AHA-/a(AH  +  HA)  +  /a^H, 

=  -(A^H  +  HA2)  +  AaH  +  /hA2  +  ai]^A  +  .  (3.23) 

''The  identities  (3.10)i  (3.11)i,  (3.13)i,  and  (3.14)2  were  observed  by  Guo  [8].  Various  authors 
have  observed  one  or  both  of  the  identities  (3.11)i,3,  often  under  the  assumption  that  A  €  Psym. 
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We  define  the  subspace  T(A)  of  Lin  as  follows: 

T(  A)  :=  {H  e  Lin  :  tr  H  =  tr  (AH)  =  tr  ( A^H)  0}  .  (3.24) 

In  particula.*".  every  tensor  H  €  T(A)  is  deviatoric,  It  iv  easily  verified  that  T(A)  = 

T(A).  By  the  Cayley- Hamilton  theorem,  H  €  T{A)  iff  tr(A*'H)  =  0  for  every 
nonnegative  integer  k\  if  A  is  nonsingular  then  H  €  T(A)  iff  tr  (A'^H)  =  0  for  every 
integer  k.  For  my  tensor  L,  let  P(L)  denote  the  subspace  of  Lin  consisting  of  all 
polynomials  in  L.  Then  T(A)  =  'P(A'*’)-'-,  the  orthogonal  complement  of  the  subspace 
of  all  polynomials  in  A”^.  Since  tr(SH)  =  0  for  every  symmetric  tensor  S  and  skew 
tensor  H,  it  follows  that 

Skw  c  T( A)  =  V{A)^  ,  V  A  €  Sym .  (3.25) 

More  generally,  suppose  that  there  is  a  basis  {ej,  62, 03}  of  V  consisting  of  eigenvectors 
of  A:  Acj  =  a,e,  (t  =  1,2,3).  Equivalently,  A  =  a.e,  (g)  e*,  where  {e*}  is  the 
reciprocal  basis  of  {Cj}.  Let  H'j  =  e*  •  He^  denote  the  components  of  the  tensor  H 

relative  to  {ej.  If  the  eigenvalues  of  A  are  distinct,  then  H  €  ^(A)  iff  =  H‘^2  — 

H^3  =  0.  If  flj  flj  =  <13,  then  H  €  T(A)  iff  =  0  and  H'^2  +  =  0-  Finally, 

=  32  =  03  =;  Cl  iff  A  =  al  iff  T(  A)  is  the  set  of  all  deviatoric  tensors. 

When  H  €  T(A),  the  identities  (3.16)-(3.18)  eind  (3.23)2  simplify  substmtially 
and  can  be  used  to  obtain  other  useful  identities.  The  results  are  summarized  in 

Proposition  3.1  Tl  e  following  identities  hold  for  any  tensor  A  and  any 
H  6  r(A)  : 


A^HA^  = 

//a AHA  -  ///a(AH  +  HA) , 

(3.26) 

A^HA  +  AHA*  = 

/aAHA  -  IlljJi , 

(3.27) 

A^H  +  HA^  +  AHA  = 

/a(AH  +  HA)  ~  //aH  , 

(3.28) 

(adj  A  )H(adj  A)  = 

7//a(/aH  -  ah  -  HA) , 

(3.29) 

AHA  = 

AaH-  (A2H  +  HA2), 

(3.30) 

A(AH  +  HA)A  = 

IIIj,H , 

(3.31) 

(adjA)H(adjA)  = 

//4(AH  +  HA). 

(3.32) 

If  A  is  nonsingular  then  (S.Sl)  and  (S.32)  are  equivalent.  If  A  is  nonsingular  then 
(3.29)  is  equivalent  to  the  identity 

///aA-^HA  *  =  /aH  -  ah  -  ha  .  (3.33) 

If  A  is  symmetric  then  (3.26)-(3.32)  hold  for  any  skew  tensor  H.  If  A  is  symmetric 
and  nonstngular  then  (3.33)  holds  for  any  skew  tensor  H. 
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Proof:  (3.26)-(3.28)  and  (3.30)  follow  immediately  from  (3.16)-(3.23)  and  the  defi¬ 
nition  of  T(A).  Then  (3.8)  and  (3.26)-(3.28)  yield  (3.29)  and  (3.33).  To  prove  (3.31), 
replace  H  with  AH  +  HA  in  (3.23)i  and  use  (3.27),  (3.28),  and  (3.11)j.  Similarly, 
(3.32)  follows  from  (3.14)i,  (3.2o),  (3.28),  and  (3.11),.  Since  T(A)  =  r(A),  (3.32) 
can  also  be  obtained  by  replacing  A  with  A  in  (3.29)  and  using  (3.9)  and  (3.10),.  If 
A  is  nonsingular  then  adj  A  =  in  which  case  (3.31)  and  (3,32)  are  easily 

seen  to  be  equivalent.  The  statements  for  symmetric  A  follow  from  these  results  and 
(3.25).  □ 

The  identity  (3.31)  is  the  key  to  our  construction  of  simple  direct  formulas  for 
Lx[H]  when  H  G  T(A).^^  Note  that  since  tr  (A^(A'"HA"))  =  tr  (A*^+'"+"H),  for 
any  nonnegative  integers  m  and  n  we  have 

A^HA"  -  A"HA'"  G  r(A),  VA,HGLin;  (3.34) 

and  if  P  G  T(A)  then  any  linear  combination  of  terms  of  the  form  A"*HA"  belongs 
to  T(A).  In  particular,  for  H  G  T(A)  each  of  the  expressions  in  Proposition  3.1 
belongs  to  T (A).  We  conclude  this  section  with  a  result  which  wdll  be  utilized  in  the 
proof  of  Theorem  4.1. 

Proposition  3.2  AHA  =  0  for  each  H  G  T(A)  iff  A  =  u  ®  V  for  some  vectors 
u  and  V.  Similarly,  AHA  =  0  for  each  H  G  T(A)  iff  A  ~  ^(u  •  v)I  —  u  (gi  v  for 
some  vectors  u  and  v.  In  both  cases,  T(A)  is  the  set  of  all  tensors  H  such  that 
tr  H  =  V  •  Hu  =  0. 

Proof:  The  “if”  part  of  the  first  result  is  straughtforward.  Conversely,  suppose  that 

AHA  =  0,  VHgT(A).  (*) 

Then  A  is  singular.  Now  in  general,  A  has  rank  zero  or  one  iff  A  =  u  g!  v  for  some 
vectors  u  and  v.  Hence,  it  suffices  to  shew  that  if  A  hais  rank  two  then  (*)  leads  to  a 
contradiction.  Since  A(b®c)  —  (b® c)A  G  T{A)  for  any  vectors  b  and  c  (cf.  (3.34)), 
(*)  implies  that 

A^b  0  A^'c  =  Ab  ®  ( AT)2c ,  Vb,cGV.  (t) 

If  A  has  rank  two,  we  may  chewse  c  so  that  A^c  ^  0  and  (A’^)2c  ^  0.  Then  (I) 
implies  that  for  every  vector  b  there  is  an  a  f  0  such  that  A(Ab)  =  A^b  =  oAb. 
Then  a  is  an  eigenvalue  of  A,  and  it  is  not  hard  to  show  that  a  is  independent  of  b, 
so  that  A2  =  oA.  If  P  :=  ^A  then  P*  =  P.  Hence,  A  =  aP  for  some  projection 

^^The  identities  (3.30)-(3.32)  were  derived  in  Scheidler  [23]  by  the  method.^  used  here  but  under 
the  assumption  that  A  is  symmetric.  A  major  special  case  of  (3.31)  was  obtained  independently 
and  by  a  different  method  by  Chen  and  Wheeler  [20].  For  symmetric  A,  they  showed  that  (3.31) 
holds  for  any  H  such  that  c'He  =  0  for  every  eigenvector  of  A.  The  set  of  all  such  H  coincides 
with  T(A)  when  A  has  three  distinct  eigenvalues  but  otherwise  is  properly  included  in  T(A). 
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P  of  rank  two.  Then  there  is  a  basis  {e^}  of  V  with  reciprocal  bajis  {e'}  such  that 
P  =  62  0  e*  -|-  63  (g)  e^,  in  which  case  (*)  implies  that  H'^2  —  ^^3  —  ^^3  ~  ^^^2  —  0 
for  each  H  e  T(A).  But  from  the  comments  preceding  Proposition  3.1,  we  see  that 
H  €  T(A)  if  ^>1  =  0  amd  //^2  =  ~H^3  /  0)  which  is  a  contradiction.  Finally,  note 
that  since  T(A)  =  T{A)  for  any  tensor  A,  we  have  AHA  =  0  for  each  H  G  T(A)  iff 
AHA  =  0  for  each  H  €  T(A)  iff  A  =  u  0  v  iff  (cf.  (3.10))  A  =  j(u-v)I  —  u0v.  □ 

4  Existence  and  uniqueness  of  solutions. 

Direct  solutions  in  T(A) 

Here,  as  in  the  previous  section,  we  do  not  require  that  A  G  Lin*,  i.e.,  that  the 
equation  AX  +  XA  =  H  has  a  unique  solution  X  for  every  H  G  Lin.  Instead,  in 
Theorem  4.1  and  Proposition  4.2  we  determine  necessary  and  sufficient  conditions  for 
the  existence  of  a  unique  solution  as  well  as  simple  direct  formulas  for  this  solution 
when  X  and  H  are  restricted  to  the  subspace  T(A).  These  results  are  then  used  in  the 
proof  of  Theorem  4.2,  which  gives  necessary  and  sufficient  conditions  for  A  G  Lin*. 
We  begin  with  the  following  result  which  is  utilized  in  the  proofs. 

Proposition  4.1  //X  G  T{A)  then  AX  -|-  XA  G  T{A);  t.e.  maps  T(A)  :nto 
itself.  Conversely,  if  A  is  nonsingular  and  AX  -|-  XA  G  T{A),  then  X  G  T(A).  In 
particular,  kerB^  C  T(A)  if  A  is  nonsi'igular. 

Proof:  The  first  part  is  ji\st  a  special  case  of  the  result  stated  after  (3.34).  Con¬ 
versely,  suppose  that  A  is  nonsingular  and  AX  -f  XA  G  T{A).  Then 
2tr(A*+^X)  =  tr(A*(AX  -I- XA))  =  0  for  any  integer  k,  so  that  X  G  T(A).'^ 
Finally,  X  G  ker  Ba  iff  AX  -|-  XA  =  0  G  ?'( A).  □ 

Theorem  4.1  The  following  conditions  are  equivalent: 

(1)  7/4^0; 

(2)  The  restriction  o/Ba  to  the  subspace  T{A)  is  nonsingular; 

(3)  For  each  H  G  '7’(A),  the  equation  AX  +  XA  =  H  has  exactly  one  solution 

X  G  T{A)  (the  possibility  of  other  solutions  outside  of  the  subspace  T (A)  is 
not  excluded). 

When  /  0,  direct  formulas  for  the  solution  X  G  T{A)  are 

lIIj,X  =  AHA  =  AHA  -  /a(AH  HA)  +  IjfH 

=  AAH-(Am-|-HA*).  (4.1) 

If  A  is  symmetric  with  IIIj^  ^  0  and  t/H  is  skew,  then  (4-1)  Is  the  only  skew  solution 
o/AX  +  XA  =  H. 

we  dropped  the  assumption  that  A  is  nonsingular  then  we  could  only  conclude  that 
tr(A*'*‘*X)  =  0  holds  for  nonnegative  k,  which  does  not  imply  trX  =  0. 
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Proof:  From  the  first  part  of  Proposition  4.1,  we  see  that  (2)<:^(3).  From  (3.31) 
with  H  — ♦  X,  we  have 

AB>,[X]A  =  A(AX-t-XA)A  =  //4X,  VXGr(A).  (*) 

If  ///^  ^  0,  then  by  (*)  it  follows  that  the  conditions  X  €  'T{X)  and  B^fX]  =  0  imply 
X  --  0.  Hence,  (l)=i>(2)'0(3);  and  (3),  together  with  (>•<),  implies  AHA  =  ///^X,  i.e., 
(4.1)i.  Then  (4. 1)2,3  follow  from  (3.23)i  and  (3.30).  If  A  is  symmetric  and  H  is  skew 
then  AHA  is  skew.  Since  Skw  C  T(A),  and  since  (4.1)  is  the  only  solution  in  T(A), 
it  follows  that  (4.1)  is  the  only  solution  in  Skw.  It  remains  to  show  that  (3)  implies 
(1).  Suppose  that  (3)  holds  and  that  =  0.  Then  by  (*),  it  follows  that  AHA  =  0 
for  each  H  G  T(A).  Hence,  by  Proposition  3.2,  we  must  have  A  =  |(u  •  v)I  —  u  (g)  v 
for  some  vectors  u  and  v,  in  which  case  H  G  T(A)  iff  tr  H  =  v  •  Hu  =  0.  But 
then  for  any  vector  w  orthogonal  to  u,  the  tensor  H  :=  u  ®  w  belongs  to  T(A)  and 
satisfies  AH  +  HA  =  0,  which  contradicts  (3).  Hence  ///^  must  be  nonzero.  □ 
Recalling  the  definition  (3.15)  of  A^,  we  see  that  (4.1)3  is  the  direct  solution  (1.16) 
obtained  by  SidorofF  [6]  and  Guo  [8]  under  the  assumptions  A  G  Psym  and  H  G  Skw. 
Various  necessary  and  suflScient  conditions  for  Illf^  ^  0  follow  from  the  results  in  the 
previous  section  and  are  summarized  in 

Proposition  4.2  The  following  conditions  are  equivalent: 

(1)  ///j^  ^0;  (2)  A  is  nonsingular; 

{3}  aU; 

(5)  lx  is  not  an  eigenvalue  of  A]  (6)  f  ^  0,  'ii  j  £  {1,2,3}. 

Since  any  nonrctil  characteristic  roots  of  A  occur  in  a  complex  conjugate  pair,  from 

(6)  we  see  that  if  A  has  a  characteristic  root  with  nonzero  real  and  imaginary  parts 
then  ^  0.  Also  note  that  A  may  be  nonsingular  even  if  A  is  singular.  Indeed, 
from  (6)  we  see  that  the  conditions 

«i=0,  03/0,  02  +  03  9^0  (4.2) 

are  sufficient  for  IIIj^  ^  0,  and  that  if  Oj  =  0  then  the  other  conditions  in  (4.2)  are 

necessary  for  ^  0.  It  follows  that  if  A  is  nonsingular  then  the  null  space  of  A 

has  dimension  at  most  one. 

By  combining  the  above  results,  we  obtain 

Theorem  4.2  The  following  conditions  arc  equivalent: 

(1)  A  G  Lin* ;  (2)  A  and  A  are  nonsingular; 

(3)  IIJx  ^  0  and  IIIj,  ^  0;  (4)  /^//^  ^  IIIx  ^0; 

(3)  IIlx  ^  0  and  IxU~^  1 1  (6)  a.  +  a,  0,  V*,;  G  {1.2,3} ; 

(7)  neither  0  nor  7^  is  an  eigenvalue  of  A. 
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Proof:  The  equivalence  of  (2)-(7)  follows  from  Proposition  4.2  and  (3.12)2.  That  (1) 
implies  (2)  follows  from  Proposition  2.1  and  Theorem  4.1.  Conversely,  if  (2)  holds  then 
by  Proposition  4.1  and  Theorem  4.1  we  have  kerB;^  C  T(A)  and  ker  B;^nT(A)  =  {0}, 
respectively.  Hence,  kerB^  =  {0}.  □ 

From  (4)  of  Theorem  4.2  we  see  that  A  €  Lin*  if  A  is  nonsinguiar  and  deviatoric. 
From  (6)  we  see  that  A  €  Lin*  if  the  characteristic  roots  of  A  have  positive  real  parts 
(e.g.,  A  €  Psym),  or  if  the  characteristic  roots  of  A  have  negative  real  parts  (e.g., 
— A  G  Psym),  or  if  A  has  a  nonzero  eigenvalue  and  a  characteristic  root  with  nonzero 
read  and  imaginairy  parts. 

The  only  part  of  Theorem  4.2  that  cao’ries  over  (with  obvious  modifications)  to 
arbitrary  dimensions  is  the  equivalence  of  (1)  and  (6).  This  result  is  a  speciad  case  of 

Theorem  4.3  Let  V  be  a  real  or  complex  inner  product  space  of  dimension  N. 
Let  A  B  €  Lin  have  characteristic  roots  a ^  and  6,  {i  =  Then  the  tensor 

equation  BX  -f  XA  =  H  has  a  unique  solution  X  for  any  given  H  iff  6,  +  ^  0  for 

each  i,j  G  {1, . . . ,  A}. 

When  V  is  a  read  TV-dimensional  inner  product  space  and  B  =  A  or  A"^,  the  condition 
6,  -f  a_,  ^  0  reduces  to  a^  -f  aj  ^  0.  As  in  the  three-dimensionad  caise,  this  condition 
cam  also  be  expressed  in  terms  of  the  TV  principad  invariants  of  A,  which  we  denote 
by  =  tr  A,  1  =  <iet  A.  Let  =  1,  let  /Jf  ^  =  0  if  ib  <  0  or  it  >  TV, 

and  let  denote  the  TV  x  TV  matrix  whose  element  in  the  ith  row  and  jth  column 
is  (-l)J'/^‘“'^\  The  matrix  (or  its  trainspose)  is  known  as  the  Hurwitz  matrix 
aissociated  with  A.  Observe  that  det  A-^  is  a  polynomiad  in  the  principad  invariamts 
of  A;  in  particular,  det  A-^  =  ~  when  TV  —  3.  For  airbitrairy  TV, 

a,  -f-  Oj  ^  0  for  each  i,j  G  {1, . .  • ,  TV}  iff  det  A^^  ^  0.  This  follows  immediately  from 
the  identities 


det  An  =  (- l)^(^+')/2  2-^'n(a, -t-a_,) 

'<> 

=  (_l)/v(^+i)/2(det  A)  Ilia,  -f  a,) ;  (4.3) 

•<j 

cf.  Halm  [25,  §2. 6,2. 7]  and  Gamtmacher  [26,  §5. 6, 5. 7]. 


*^Most  of  the  proofs  in  the  literature  deal  only  with  the  special  case  where  V  is  the  vector  space 
of  N-tuples  of  complex  numbers  and  A,  B,  H,  and  X  are  complex  N  x  N  matrices;  cf.  Gantmacher 
[29,  §8.3],  Bellman  [31,  §12.13],  Jacob  and  Polak  [32],  and  Feintuch  and  Rubin  [33]  for  four  different 
proofs.  These  proofs  utilize  the  fact  that  every  characteristic  root  of  A  is  an  eigenvalue  of  A;  hence, 
they  do  not  CArry  over  directly  to  a  real  vector  space.  However,  the  corresponding  theorem  for  real 
matrix  equations  and,  consequently,  the  general  theorem  stated  above,  can  be  obtained  from  the 
result  for  complex  matrix  equations. 
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5  Direct  formulas  for  L^,  and 
in  three  dimensions 

In  this  section  we  assume  thci  A  G  Lira*  and  V  is  three-dimensional.  Then  from 
Theorem  4.1  and  (1.18)  it  follesvs  that  for  each  H  G  T(A), 

IllxUiil]  =  AHA 

-  AHA  - /a(AH  +  HA)  + 

=  AaH  -  (A2H  +  HA2) ;  (5.1) 

in  particular,  these  formulas  hold  whenever  ±A  G  Psym  and  H  G  Skw,  in  which  case 
La[H]  is  skew.  For  arbivraxy  H,  all  of  the  direct  formulcis  for  La[H],  Ma[H],  aad 
Na[H]  derived  below  will  be  obtained  from  (5.1)i_2,  the  results  in  Sections  2  and  3, 
and  the  fact  that  AH  —  HA  G  T(A)  for  each  A,  H  G  Lin,  which  is  just  a  special  case 
of  (3.34).  This  last  fact  allows  us  to  replace  H  with  AH  —  HA  in  (5.1). 

By  replacing  H  with  AH  -  HA  in  (5.1)i  and  using  =  La[AH  —  HA]  cind 

A  =  IjJ.  —  A,  we  obtain 

///aMa[H]  A(AH  -  HA)A  =  A(HA  -  AH)A 

=  A(AHA)  -  (AHA)A  =  (AHA)A  -  A(AHA1 
=  AHA2-A2HA  (5.2) 

for  ciny  tensor  H.  Similarly,  from  (5.1)2  with  H  — »  AH  —  HA,  we  obtain 

///aMa[H]  =  A^H  A  -  AHA2  +  /^(HA^  -  A^H)  +  /a*(AH  -  HA) .  (5.3) 

A  lengthier  derivation  of  this  result  follows  by  replacing  H  with  AH  —  HA  in  (5.1)3; 
this  yields  some  A^  terms  which  can  be  reduced  by  means  of  the  Cayley-Hamilton 
theorem  (3.7).  Mehrabadi  and  Nemat-Nasser  [11]  obtained  (5.3)  by  repeated  appli¬ 
cations  of  the  Cayley-Hamilton  theorem. 

In  view  of  (1.19),  the  direct  formulas  (5.2)  and  (5.3)  yield  a  variety  of  direct 
solutions  X  =  Ma(H]  of  the  tensor  equation  AX-f  XA  =  AH  — HA,  H  arbitrary.  In 
particular,  if  A  and  H  are  symmetric  then  Ma(H]  is  skew,  and  we  have  the  cdternate 
formulas 


and 


7//aMa[H] 


2  skw  (A  AH  A)  =  2skw(AHA2) 

2  skw  (a^HA  -  /aA^H  -f  /a^AH) 


-///aMaI'h; 


=  2  skw  (AHA A)  =  2  skw  ( A^H A) 

=  2 skw  (aHA^  -  /aHA2  +  IaHa)  . 


(5.4) 


(5.5) 


20 


Next,  we  derive  direct  formula.s  for  Lx[H]  and  Nj^[H]  for  arbitrary  H.  From 
(2.18)  and  (5.2)i,  we  obtain 


2U[H]  = 


H  + -^A(HA  -  AH)A 
1 


=  HA-‘ +  — A(H- AHA-i)A 
=  A-»  H  + -^A(AH  -  HA)A 

-  A-»H  +  -^A(H- A-iHA)A. 
■“■'A 


(5.6) 


Other  formulais  for  Lj^[H]  follow  from  (2.18),  (5. 2)2-5,  and  (5.3)-(5.5).  In  particular, 
if  A  and  H  are  symmetric  then  L^fH]  is  symmetric,  and 


La[H]  = 


+  ^skw  (A^HA; j  A  1 


=  A  -1  I  +  7^-skw  (AHA2) 
2  A“A 


(5.7) 


Alternate  expressions  for  L^fH]  follow  by  substituting  (8.8)2  for  A-r. 

Since  Na[H]  =  M^fAH  -  HA)  =  Ma[HA  -  AH],  by  replacing  H  with  AH  -  HA 
in  (5.2)1  and  replacing  H  with  HA  —  AH  in  (5.2)2,  ''''^0  obtain 

=  A(A2H  +  HA2  -  2AHA)A 
=  A2(AHA)  +  (AHA)A2  -  2A(AHA)A 
=  A(A2H  +  HA2 -2AHA)A 

=  A^HA  +  AHA^  -  2A2HA2 .  (5.8) 


Another  expression  for  N^fH]  can  be  obtaiined  by  using  A  =  l^l  —  A  in  (5. 8)1, 
expanding,  and  then  using  (3.7)  to  reduce  the  A^  terms.  Similarly,  we  can  reduce  the 
A'’  terms  in  (5. 8)4  by  (3.7)  with  A  A.  The  results  are 


=  -2A2HA2  +  2/a(A2HA  +  AHA^)  -  2/4AHA 
+  Fa  (AH  +  HA)  - 

=  -2.A2HA2  +  /^(A^HA  +  AHA^)  -  2//aAHA 

+  ///a(AH  +  HA) ,  (5.9) 

where  alternate  expressions  for  are  given  by  (3.13),  and 

r*  = /aWa  +  w/a  = ///a  +  2/w.. .  (5.10) 
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By  using  Rivliu’s  identities  (3.16)  and  (3.17)  in  (5.9)i,  we  obtain 

=  -4//aAHA  +  [Ik^Ik  +  3///a)(AH  +  HA)  -  4/a///aH 

-  2  (a™„  A=  +  A  +  ///*£.S;;„l)  ,  (5.  n ) 

where  ^  and  cta^h  given  by  (3.20)  and  (3.21),  and 

»7a,h  “  “  -^A^A^H  +  Aa/ah  +  ■  (5.12) 

Finally,  by  using  Rivlin’s  identity  (3.18)  in  (5.11),  we  obtcxin 

///aNa[H]  =  4//A(A2H  +  HA2)-3//4(AH  +  HA)-f  4(/7a"-/a^//a)H 

“  2  (va.h^^  +  '/a.rA  +  »?a,hI)  >  (5.13) 

where 

^A,H  =  —  -^A^AH  +  2//a-^H  >  (^''•14) 

r}!]^  =  +  (.5.15) 


=  2//a/a^h  +  (///a  -  2/a//a)/ah  +  (2//a^  - /a///a)  4  •  (5-16) 

In  the  remainder  of  this  section  we  derive  some  additional  formulcis  for  La[H] 
for  arbitrary  H.  By  substituting  the  formula  (5.9)i  for  Na[H]  into  (2.19)i  and  then 
using  (3.23)i,  we  obtain 

2//4La[H]  =  AHA  - /^(AH  +  HA)  + /4H 

+  /a///aA-^HA-^  -  7//a(A-JH  +  HA-i) 

=  AHA  +  77aH 

+  7a777aA-’HA-»  -  777a(A-‘H  -b  HA'^)  •  (5.17) 

The  direct  formula  (5.17)i  (with  777a  replaced  by  the  equivalent  expres¬ 

sions  (3.11)1  and  (3.13)i)  was  stated  without  proof  by  Leonov  [5]  cind  Stickforth 
and  Wegener  [12]  for  the  case  A  €  Psym  and  H  €  Sym.^*  If  we  use  the  relation 
777aA~^  =  adj  A  in  (5.17),  we  obtain 

27/7ji^LA[H]  =  AHA  -  7a(AH  +  HA)  +  77aH 

-t-  ^  (adj  A)H(adj  A)  —  (adj  A)H  —  H(adj  A) 

777a 

=  AHA-t-77AH 

+  7~(adj  A)H(adj  A)  -  (adj  A)H  -  H(adj  A) .  (5.18) 

*'7iA 

‘“^Leonov  attributes  the  result  to  L.  M.  Zubov.  Stickforth  and  Wegener  refer  the  reader  to  some 
lecture  notes  by  Stickforth. 
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If  the  formula  (8.8)2  substituted  into  the  last  term  in  (5.17)i,  we  obtain 

2///aLa[H]  =  AHA  -  (A^H  +  HA*)  -f  AaH  +  /a///aA-*HA-i  (5  19^ 

The  formulas  (5.18)2  and  (5.19)  are  due  to  Muller  [16]  and  Jameson  [15],  respectively.'® 
By  using  (3.8)  in  (5.19)  or  (5.18)i  and  expanding,  we  obtain 

2//.rA///ALA[H]  =  /aA*HA*  -  /a*  (A*HA  +  AHA*)  +  HIj,  (A*H  +  HA*) 

+  +  JIU)  AHA  -  /a^//a(AH  +  HA)  +  qaH  ,  (5.20) 

wh.crG 

Oa  =  /a'/Z/a  + //a/Wa  -fA^A  +  W/aAa 

=  ZaZZa  +  ZaVZZa  -  ZZaZZZa  .  (5.21) 

In  view  of  (1.18)  and  (3.11)i,  this  is  the  direct  formula  (1.17)  of  Hoger  and  Carlson 
19];'^'  cf.  also  Smith  [14].'® 

Substitution  of  Rivlin’s  identities  (3.16)-(3.18)  into  (5.20),  together  v/ith  (3.23)i, 
yields 


///aLa[H]  =  AHA  - /a(AH  +  HA)  + /a^H 
+  +  A‘,’hA  + 

=  AHA  +  A^  +  A',h  A  +  A°.hI  .  (5.22) 

where 

2/^A,H  =  JJ^  ~  ^A^AH  +  ^A-^H  ~ 

=  jjJ2  (“^a-^Aah  +  =  -^Aa-ih  ’  (5.23) 

'®*4uller  and  Jameson  obtained  direct  solutions  of  the  tensor  equation  BX+XA  =  H  for  arbitrary 
dimensions.  Muller’s  formula  is  quite  cor^iplicated.  For  the  raise  B  =  A''^,  both  authors  simplified 
their  general  results,  and  Jameson  listed  the  corresponding  formulas  in  two  and  three  dimensions. 
For  B  =  A^  and  dim  V  =  3,  the  formulas  of  Muller  and  Jameson  are  equivalent  to  those  obtained 
irom  (5.18)2  and  (5.19),  respectively,  via  Proposition  2.4;  they  redur  to  (f  ’8)2  and  (5.19)  when 
A  is  symmetric.  Their  general  formulas  (for  A  and  P  unrela'  d)  also  reduce  to  (5.18)2  and  (5.19) 
when  B  =  A. 

'^Hoger  and  Carlson  showed  that  (1.17)  is  a  solution  of  AX  +  XA  =  H  for  any  tensor  H  and  any 
tensor  A  such  that  IIIa{IaOa  -  HIa)  ^  0.  They  established  uniqueness  of  this  solution  when  A  is 
also  symmetric. 

'*For  symmetric  H,  Smith  obtained  a  direct  tiolution  of  the  tensor  equation  A^X  +  XA  =  H 
for  arbitrary  N  —  dimV.  This  solution  has  the  form  (det/l7<)X  =  X)mn=i  Cm,n(A^  )'"“'HA^~', 
where  Ati  is  the  Hurwitz  matrix  associated  with  A  (cf.  Section  4),  and  the  coefficients  c^.n  are 
complicated  polynomials  in  the  principal  invariants  of  A  involving  cofactors  of  A-h.  Smith  gave 
explicit  expressions  for  these  coefficients  when  di.iiV  =  2,3,4.  For  dimV  =  3,  his  formula  is 
equivalent  to  the  one  obtained  frorri  (1.17)  via  Proposition  2.4;  this  formula  reduces  to  (1.17)  when 
A  is  symmetric. 
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=  JJj-  -^A  +  {^A  -  I^^a)  ^AA  +  ^A^A 

—  ~^AH  +  2/x/h  ~  Ia^a-^h  —  ^Ah  ”  ^a^Aa-^h  '  (5-24) 

=  JjJ-  +  0A^AH  +  Pa^h) 

=  -^a^Ah  +  ^^^A^A-‘H  "=  ^^^a(^a^a-’h  “  -^a->h)  ^  (5.25) 

0A  =  /a(///a-///a)  =  -^a(2///a-/a//a),  (5.26) 

=  IUIIIa-  IaIIIa  =  IaIIa  -  I^UIIa 

=  IaIIa  -  IaHU  -  IIaIIIa  •  (5-27) 

From  (5.22)  and  (3.23)2,  we  also  have 

///aLx[H]  =  -(A2H  +  HA2)  + AaH 

+  7^?hA^  +  7i\^HA  +  ,  (5.28) 

where 


27a, H  =  Jjj^  i^A^Am  ~  Ia^ah  +  r a.^h) 


-  (  Ia^AaH  +  Ta^h)  -  .^H  d-  Ia^A-^A  1 

(5.29) 

^7a,h 

-  [  ^A I  Am  +  (.^A  +  in  a)  Iaa  I  a  IIaIa] 

=  .^AH  ~  lA^lA-m  =  —{IaA  +  lAlAA-m)  ’ 

(5.30) 

27a,H 

-  jjj  {l^AlAm  IaIIaIaa  +  KaIa) 

-■  -IaIaa  +  Ta^a-ih  =  77/a(/a-ih  +  ^a^a-»h)  ■> 

(5.31) 

Fa  is  given  by  (5.10),  and 

7a  =  IIaTa-  IaIIIa  =  IaIJa!  -  IIIa^a 

=  /a//a^-/a"///a  +  //a///a-  (5.32) 

The  direct  formula  (5.28)  with  the  expressions  (5.29)3,  (5.30)2,  and  (5.31)3  for  7a(h  i® 
due  to  Sidoroff  [6].*'’  The  direct  formula  (5.28)  with  the  expressions  (5.29)i,  (5.30)i, 
and  (5.31  )x  for  7a|h  was  obtained  by  Roger  and  Carlson  [9]  by  essentially  the  same 
method  as  used  here,  i.e.,  from  the  direct  formula  (5.20)  and  Rivlin’s  identities.  Roger 

^®For  A  €  Psym,  Sidoroff  derived  (5.1)3  for  skew  H  and  (5.28)  for  symmetric  H  and  then  conjbined 
these  to  obtain  (5.28)  for  aibitreuy  H.  He  noted  that  (5.28)  yields  a  solution  of  AX  +  XA  =  H  for 
any  tensor  A  such  that  IIIa{IaOa  -  HIa)  0>  but  did  not  establish  uniqueness  for  this  case. 


2A”>h 


and 
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and  Carlson  observed  that  (5.28)  collapses  to  the  formula  (5.1)3  of  Sidoroff  and  Guo 
when  A  €  Fsym  and  H  is  skew*.  Indeed,  for  any  A  €  Lin*  and  any  H  G  T{A),  we  see 
that  (5.22)  and  (5.28)  collapse  to  (5.1),  2  and  (5.1)3,  respectively.  The  direct  formulas 
(5.6),  (5.7),  and  (5.17)-(5.20)  do  not  have  this  property. 

From  (5.22)2  and  the  identity  A*  =  A^  +  2/xA  —  1x1,  we  obtain 

OTaUIH]  =  A(m  AV)A-71','hA  +  <1;,’hI 

=  A(H  +  AV)A  +  7i',',iA  +  4V.  (5-33) 

where 

24"’h  “  =  (5.34) 

25a, H  ==  {5a^a»h  +  5a3ah  +  ^^a5a.^h) 

=  1x1a.h  +  >  (5.35) 

and 

iA  =  nil,  -  /a  =  -(/aA.a  +  IIIk)  =  -k  (2^a  +  /a-)  ,  (5.36) 

4  =  /a  +  4  =  .fA(fflA  -  5a)  - /a(/aAa  +  2W/a) 

=  /a(/a  -  IkIU  +  2///a)  .  (5.37) 

Also,  by  (5.23)  and  (5.30),  we  have 

■“  Ta.H  ~  -^A^A.H  +  2-^Ah  '  (5.38) 

We  conclude  this  section  by  considering  the  special  case  where  A  and  H  are 
time-dependent  tensors  with  H  related  to  A,  its  time  rate  of  change  A,  and  some 
time-dependent  tensor  Z  as  follows: 

H  =  A  4  AZ  -  ZA ;  (5.39) 

cf.  (6.4)  and  (7.14).  Then  tr(p(A)H)  =  tr(p(A)A),  where  p(A)  is  any  polynomial 
in  A.  In  particular,  we  have 

/h  =  /a  =  4  ,  /ah  =  //a  ,  (5.40) 

/ah  =  /aA  =  —  /a /a  ~  //a  »  (5.41) 

/a*h  =  /a*A  =  =  /a /a  -  ///a  ,  (5.42) 

/(adj  A)H  =  /(.djA)A  =  10 X ,  /a->h  =  /a- A  =  ///a////a  .  (5.43) 

These  identities  yield  alternate  expressions  for  the  coefficients  in  the  direct  formulas 
(5.22),  (5.28),  and  (5.33)  for  L2,^[H]  and  for  the  coefficients  in  the  identities  (3.16)-- 
(3.18)  and  (3.23)2. 
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6  Applications  to  kinematics  of  continua 

In  this  section  we  apply  the  preceding  results  to  the  derivation  of  some  kinematic 
formulas  for  a  material  undergoing  a  smooth  motion  in  three-dimensional  space.  By 
taking  the  material  time  derivative  of  the  relations  =  C  and  =  B,  we  obtciin 
the  tensor  equations  (1.6)  for  U  and  V.  Hence, 

U  =  Lu[C]  and  V  =  Lv[B].  (6.1) 

A  variety  of  direct  formulas  for  U  in  terms  of  C  and  U,  and  for  V  in  terms  of  B  and 
V,  follow  from  (6.1)  and  the  direct  formulas  for  Lx[H]  in  Section  5.  In  particular, 
the  formulas  that  follow  from  (5.20)  and  (5.28)  are  due  to  Hoger  and  Carlson  [9]. 

By  differentiating  B  =  FF"^  and  using  the  relation  between  the  velocity  gradient 
and  the  material  time  derivative  of  the  deformation  gradient, 

F  =  LF,  (6.2) 

we  obtain  the  following  formulas  for  the  material  time  derivatives  of  B  and  B“M 

B  =  LB-hBLT  and  -  (B-» )  =  B-‘BB-' =  B-»L  +  L^B-^  (6.3) 

Then  (6.3),  (1.4),  and  the  definition  (1.5)  of  the  Jaumann  rate  yield  the  tensor  equa¬ 
tions  (1.10)  for  D.  Hence, 

D  =  Lb[6)  =  -Lb-.1(B->)“1.  (6.4) 

A  variety  of  direct  formulas  for  D  in  terms  of  6  and  B,  or  in  terms  of  (B“^ )°  and  B~* , 
follow  from  (6.4)  and  the  direct  formulas  for  L^[H]  in  Section  5.  In  particulair,  the 
formulas  which  follow  from  (6.4)  and  (5.17)i  are  due  to  Leonov  [5],  and  the  formula 
which  follows  from  (6.4)}  and  (5.28)  is  due  to  Sidoroff  [6].  Also  note  that  the  identities 
(5.40)-(5.43)  can  be  used  in  the  expressions  for  the  coefficients  in  the  direct  formulas 
for  D  which  follow  from  (6.4),  (5.22),  (5.28),  and  (5.33). 

Let  Lji  denote  the  rotated  velocity  gradient'. 

La  :=  R^LR  =  Dr  +  Wa  .  (6.5) 

Here  Da  and  Wa  denote  the  rotated  stretching  and  spin  tensors: 

Da:=R'^DR  =  symLa  and  Wa  :=  R'^WR  =  skw  La .  (6.6) 

By  difft  entiating  C  =  F^F  and  using  (6.2),  (1.4)2,  F  =  RU,  and  (6.6)i,  we  obtain 

C  =  2FTDF  =  2UDbU  .  (6.7) 
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This,  together  with  (1.6)i,  yields  the  tensor  equations  (1.9)  for  U.  Hence, 

U  =  2Lu[UDaU]  =  2Lu-.[Dr]  .  (6.8) 

Although  direct  formulas  for  U  in  terms  of  Dr  and  U  follow  from  (6.8)  and  the 
direct  formulas  for  L^.  [H]  in  Section  5,  we  can  obtain  simpler  formulas  by  using  the 
identities  (2.21).  Indeed,  from  (6.8)  and  (2.21)  we  obtain  the  following  relations  for 
the  material  time  derivative  of  U: 


U  =  U(Da-M„[Da])  =  (Da  +  M„|Dtt|)U 

=  KUDb  +  DaU  -  N„lDa)) .  (6.9) 


Then  a  variety  of  direct  formulas  for  U  follow  from  (6.9)i,2  direct  formulas 

(5.2)-(5.5)  for  Mx[H].  In  particular,  we  have 


U 


Dr  + 
Dr  + 
Dr  -f 


1 


Ilh 

1 

nii 

2 


U(UDr-DrU)U  u 
(UDrU*  -U2DrU)|  u 


skw(UDaU2) 


u. 


(6.10) 


Some  of  the  formulas  which  follow  from  (6.9)3  and  the  direct  formulas  for  N^fH]  in 
Section  5  are 

2  U  =  UDr  +  DrU  +  ^U(2UDrU  ~  U^Dr  -  DrU2)U  ,  (6.11) 

lllfj 

and 

///tjU  =  U^DrU^  - /u(U2DaU  +  UDrU*)  + //^UDrU 
-  ///u(UDr  +  DrU)  +  /u///uDr 
=  27/uUDaU-2///u(UDR  +  DRU)  +  2/u///uDR 
+  +  IIIv 

=  -2//u(U2Dr  +  DrU2)  +  2///o(UDr  +  DrU) 

+  2(/u7//u  -  -^-^u  )Dr  +  »7u,d,iU^  +  »7u,DrU  -f  Vu.DrI?  (6.12) 

where  the  q6)  axe  defined  by  (3.19)-(3.22)  and  the  r?6)  axe  defined  by  (5.12)  and 
(5.14)-(5.16).  The  direct  formulas  (6.12)i  3  were  obtained  by  Mehrabadi  and  Nemat- 
Nasser  [11]. 20  The  formula  (6.12)3  ^Jso  derived  by  Guo,  Lehmann,  and  Liang 
[13].  Now  for  any  integers  m  and  n. 


U"‘  =  RTV'"R  and  U^-DrU"  =  RT(V"*DV")R  .  (6.13) 

^°Their  formula  corresponding  to  (6.12)i  contains  a  misprint;  ZZu  should  be  replaced  with  J27u 
in  the  last  term  of  their  formula. 
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By  use  of  (6.13)  we  may  convert  (6.10)-(6.12)  to  formu,a.s  for  U  in  terms  of  D,  V, 
and  R.  When  this  conversion  is  applied  to  (6.12),,  we  recover  a  formula  obtained  by 
Roger  [10]. 

From  (1.6)2,  (6.3)i,  and  B  =  V^,  we  obtain  the  tensor  equation  (1.7)  for  V. 
Hence, 

V  =  Lv[V2LT  +  LV2].  (6.14) 

Direct  formulas  for  V  in  terms  of  L  and  V  follow  from  (6.14)  and  the  direct  formulas 
for  Lx[H]  in  Section  5;  in  particular,  the  formula  for  V  obtained  from  (5.17),  is 
due  to  Stickforth  and  Wegener  [12].  However,  simpler  formula:?  can  be  obtained  by 
using  the  identities  (2.26).  Indeed,  from  (6.14),  (2.26),  (1.4),  and  (1.5),  we  obtain  the 
following  relations  for  Jaumann  rate  of  V: 

V  =  VD  -  Mv[D]V  =  DV  +  VMv[D] 

=  1(VD  +  DV  4- Nv[D]) .  (6.15) 

These  relations  can  also  be  derived  as  follows.  From  (1.7)  and  (1.4),  we  have 

VV  +  VV  =  V^D  +  DV2  -  V^W -f  WV2 .  (6.16) 

o 

But  this  is  easily  seen  to  be  equivalent  to  the  tensor  equation  (1.8)  for  V.  Hence, 
V  =  Lv[V2D  +  DV2],  and  the  identities  (2.24)  yield  (6.15).  A  variety  of  direct 
formulas  for  the  Jaumann  rate  of  V  follow  from  (6.15),  2  and  the  direct  formulas  for 
Ma[H1  in  Section  5.  In  particular,  we  have 

/7/vC^-DV)  =  VV(VD-DV)V 
=  V(VDV2-V2DV) 

=  Vskw(VDV2).  (6.17) 


Some  of  the  formulas  which  follow  from  (6.15)3  and  the  direct  formulas  for  Njiv[H]  in 
Section  5  are 


and 


2^  =  VD  +  DV  +  -^V(V2D  +  DV2  -  2VDV)V  , 

Ill-fy 

=  -V^DV^  +  /v(V2DV  +  VDV2)  -  IlyYDV 
+  /v//v(VD  +  DV)  -  /v///vD 
=  -V2DV2  ^  /y(v2DV  +  VDV2)  -  7/^VDV  +  7v777vD 


(6.18) 


=  -277vVDV  -}-  (7v77v  +  777v)(VD  +  DV)  -  27v777vD 
-  «v!dV^  -  -  777vay;i,i 

=  277v(V2D  +  DV2)  -  777v(VD  +  DV)  +  2(77v^  -  7v777v)D 


(6.19) 
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The  formula  (G.IQ)^  was  obtained  by  Guo,  Lehmann,  and  Liang  [13]. 
Next,  we  derive  direct  formulas  for  the  skew  tensors 


0  :=  RRT  =  -RRT  =  RfijiRT 

(6.20) 

and 

Or  ;=  R1 R  =  -RTR  =  RT^R  . 

(6.21) 

Then  direct  formulas  for  R  follow  from  the  relations 

R  =  riR  =  Rf^R  . 

(6.22) 

By  differentiating  the  polar  decompositions  (1.2)  and  using  (6.2),  (6.20)i, 
the  fact  that  V  =  V'^  and  U  -  U'^,  we  obtain  the  formulas 

(6.21)1,  and 

V  -  LV  -  vn  =  VLT  +  nv 

(6.23) 

and 

U  =  (La  -  na)U  =  U(LT  +  Ha) , 

(6.24) 

where 

LT  ;=  (Lr)T  =  RTLTR  =  Dr  -  Wr  . 

(6.25) 

Observe  that  (6.23)2  (6.24)2  equivalent  to  the  following  equations  for  the  skew 

tensors  ft  and 


vn  +  nv  =  L  and  Ulla  +  =  Lr  ,  (6.26) 

where 

L  :=  LV- VLT  =  2skw(LV) 

=  DV-VD  +  WV  +  VW,  (6.27) 

and 

Lr  :=  LrU-ULT  =2skw(LaU)  =  RTLR 

=  DrU-UDr  +  WrU  +  UWr.  (6.28) 

In  particulcir,  (6.26)i  with  L  given  by  (6.27)i  is  the  tensor  equation  (1.11)  for  fl. 
From  (6.26)  it  follows  that 


=  LyfL]  and  ^r  —  L.u[Lr]  .  (6.29) 

Note  that  the  tensor  equation  (6.26)2  can  also  be  obtained  by  multiplying  the  terms 
in  (6.26)i  on  the  left  by  R^  and  on  the  right  by  R,  and  then  using  U  =  R^VR, 
(6.21),  and  (6.28).  Similarly,  (6.29)2  also  follows  from  (6.29)i  by  using  (2.15)  with 
K  =  L,  A  =  V,  and  Q  =  R'^.  Since  L  and  Lr  are  skew,  direct  formulas  for  Q  and 
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Oh,  follow  from  (6.29)  and  the  direct  formulas  (5.1)  for  Lji^[H].  In  particular,  for  O 
we  obtain 


///^O  =  VLV  =  VLV  -  Jv(VL  +  LV)  +  VL 

=  AvL-(V2L-fLV2)  =  AvL-(BL  +  LB).  (6.30) 

Various  expressions  for  the  invariant  Ay  follow  from  (3.15);  in  particular,  we  have 

Ay  :=  /v  “  -ffv  ~  Ib  -f-fv  =  ■  (6.31) 

The  direct  formula  (6.30)3  for  O  is  due  to  Guo  [8].  The  formula  (6.30)i  was  obtained 
recently  by  Chen  and  Wheeler  [20].^^ 

Prom  (6.30)  and  (6.27)  we  see  that  O  depends  on  L  euid  V,  or,  equivalently,  on 
D,  W,  and  V.  Since  W  and  O  are  both  tensor  measures  of  material  spin,  one  might 
expect  the  dependence  of  f2  on  W  to  be  of  a  simpler  form  than  its  dependence  on  D 
and  V.  This  turns  out  to  be  the  case,  although  it  is  not  obvious  from  the  formulas 
(6.30).  Instead,  we  use  the  identity  (2.27),  together  with  (6.29)i,  (6.27)i,  (1.4),  and 
(2.23),  to  obtain  the  relations 

n  =  W-My[D]  =  W  +  My-,[D].  (6.32) 

(6.32)1  may  also  be  obtained  as  follows.  Rrom  the  expression  (6.27),3  for  L,  we  see 
that  (6.26)i  is  equivalent  to  the  tensor  equation  (1.12)  for  W  — then  by  (1.19)  've 
have  W  -  n  =  My[D].  FVom  (6.32)i  and  the  direct  formulas  for  in  Section  0, 

it  follows  that 

///v(n-W)  =  V(DV- VD)V  =  2skw(VDVV) 

=  V( VD  -  DV)V  =  V^i^v  -  VDV2  =  2  skw  (V^DV) 

=  VDV2  -  V=DV  +  /y(V2D  -  DV2)  +  V(DV  -  VD) 

=  2  skw  ( VDV2  -  /yDV2  +  VDV)  .  (6.33) 

Since  W  is  skew,  we  may  also  obtain  (6.33)i  from  (6.30)i,  the  expresssion  (6.27)3 
for  L,  and  the  identity  (3.31).  From  (6.27)-(6.29)  we  see  that  results  cinaJogous  to 

had  also  derived  (6.30)i  prior  to  seeing  their  paper.  Chen  and  Wheeler  [20]  first  derived  a  direct 
formula  for  DpR,  the  derivative  of  the  rotation  tensor  with  respect  to  the  deformation  gradient, 
and  then  obtained  (6.30)i  from  the  relations  R  =  DpR[F]  and  F  =  LF.  Although  they  did  not 
solve  any  tensor  equations  in  deriving  their  formula  for  DpK,  they  did  utilize  a  special  ca£>e  of  the 
identity  (3.31);  cf.  the  footnote  in  Section  3. 

^^The  tensor  equation  (1.12)  was  stated  (without  proof)  by  Green  [34]  and  derived  by  Green  and 
Meinnb  [35].  In  both  papers  it  was  noted  that  detV  ^  0  and  stated  (without  proof)  that  this 
condition  implies  that  (1.12)  may  be  solved  for  O  —  W  in  terms  of  V  and  D,  but  no  solutions  were 
given.  The  tensor  equations  (1.7),  (i  ll),  and  (1  12)  were  derived  by  Dienes  [7],  He  obtained  the 
solution  w  =  w4-  of  (1.12),  where  m,  w  and  s  are  the  axial  vectors  of  fl,  W,  and  DV  -  VD, 
respectively,  but  did  not  obtain  an  explicit  formula  (cf.  (3.14))  for  V~^ 
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(6.30)-(6.33)  and  (1.12)  also  hold  for  They  may  be  obtained  from  (6.30)-(6.33) 
and  (1.12)  by  the  replacements 

V  U,  B  ->  C,  ^  Qr,  L  -4  Lr,  W  -4  Wr,  D  -  Dr  .  (6.34) 

Hoger  [10]  derived  (6.33)6  by  substituting  her  formula  for  U  into  the  identity 

n-W  =  RU-iURT-D,  (6.35) 

which  follows  from  (6.24)3.  Stickforth  and  Wegener  [12]  derived  a  formula  equivalent 
to  (6.33)7  by  substituting  their  formula  for  V  into  the  identity 

D  =  V-^LV- V-iy,  (6.36) 

which  follows  from  (6.23)i.  Our  derivation  of  (6.83)5  is  closer  in  spirit,  though  different 
in  detail,  to  that  of  Mehrabadi  and  Nemat-Nasser  [11].  They  derived  the  analog  of 
the  tensor  equation  (1.12)  for  Hr  and  solved  this  equation  to  obtain  the  analog  of 
(6.33)6  for  Dr;  they  also  iiCvcd  the  corresponding  result  for 

A  vaxiety  of  additi.ui  1  formulas  for  V  and  U  can  be  obtained  by  substituting  the 
formulas  (6.30)  for  fl  into  (6.23),  and  by  substituting  the  analog  of  (6.30)  for  Dr  into 
(6.24).  In  particular,  the  formulas  for  V  and  U  which  follow  from  (6.86)3,  (6.23)i, 
and  (6.24)i  are  due  to  Guo  ^8],  and  the  formula  for  V  which  follows  from  (6.30)i  and 

(6.23) i  was  obtained  by  Chen  and  Wheeler  [20].  From  (6.15)i_.3  (6.32)i,  or  from 

(6.23) ,  (1.4),  and  (1.5),  we  obtain  the  following  simple  expressions  for  the  Jaumann 
rate  of  V  in  terms  of  D,  D  —  W,  and  V: 

f  =  VD  +  (D-W)V  =  DV-V(n- W) 

=  i[VD  +  DV  +  (D- W)V- V(D -W)] 

=  sym  [(D  +  D  -  W)V] .  (6.37) 

By  substituting  the  formulas  (6.33)  into  (6.37),  we  recover  the  formulas  (6.17)  and 
(6.18).  Similarly,  from  (6.9)i_2  fbe  analog  of  (6.32)]  for  Dr,  or  directly  from 

(6.24)  and  (6.5),  we  see  that 

U  =  [Dr  -  (Dr  -  Wr)]U  =  U(Dr  +  Dr  -  Wr) 

=  i[UDR  +  DrU  +  U(Dr  -  Wr)  -  (Dr  -  Wr)U] 

=  sym[D(DR  +  OR- Wr)].  (6.38) 

The  formulas  (6.38)13  were  obtained  by  Mehrabadi  and  Nemat-Nasser  [11].  By  sub¬ 
stituting  the  analog  of  (6.33)  for  Dr  —  Wr  into  (6.38),  we  recover  the  formulas  (6.10) 
and  (6.11). 

^^Equations  similar  to  (1.12)  arise  in  the  kinematics  of  elastic-plastic  deformations  based  on  the 
multiplicative  decomposition  of  the  deformation  gieidient  into  elastic  and  plastic  parts;  cf.  Nemat- 
Nasser  [36,  37]  and  Obata  e<  al.  [38]. 
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7  Additional  kinematic  formulas 


If  the  deformation  gradient  F  is  known  then  the  Canchy-Green  tensors  C  =  F^'^F 
and  B  =  FF"^  are  easy  to  calculate,  whereas  the  stretch  tensors  U  =  \/C  and 
V  =  ^e  generally  more  difficult  to  compute  since  they  involve  a  tensor  square 
root.  Hence,  for  some  applications  it  might  be  useful  to  have  direct  formulas  similar 
to  those  in  Section  6  but  involving  C  or  B  instead  of  U  or  V.  To  this  end  we  utilize 
a  formula  for  the  square  root  of  a  symmetric,  positive-definite  tensor  due  to  Ting  [39] 
and  Stickforth  [30] 


///v V  =  - B2  -I-  A vB  -I-  /v///vl ,  (7.1) 

with  an  analogous  formula  for  U  in  terms  of  C.  Recall  that  Av  is  given  by  (6.31), 
and  note  that 

77v  =  |(V-/b),  777v=y^,  777v  =  7v77v-  777v.  (7.2) 

Thus  if  we  wish  to  express  V  solely  in  terms  of  B  and  its  principal  invariants,  we 
need  an  explicit  formula  for  ly  in  terms  of  the  principal  invariants  of  B.  In  three 
dimensions  this  formula  is  rather  complicated.  Of  course,  given  B  one  could  also 
compute  the  eigenvalues  of  B  and  then  compute  the  coefficients  in  (7.1)  from  the 
eigenvcdues  of  V.  In  any  case,  by  substituting  (7.1)  or  its  analog  for  U  into  the 
direct  formulas  in  Section  6,  we  obtain  corresponding  formulas  in  terms  of  B  or  C. 
The  scalar  coefficients  in  these  new  formulas  are  more  complicated  and  involve  the 
principal  invariants  of  IT  or  V;  cf.  Hoger  and  Carlson  [9],  where  direct  formulas  were 
obtained  for  U  in  terms  of  C  emd  C,  and  for  V  in  terms  of  B  and  B.  Here  we  list 
only  the  simplest  of  the  results  which  can  be  derived  by  this  procedure,  namely  the 
formula  for  fl  -  W  obtained  from  (7.1)  and  (6.33)6; 

777v"(fl-W)  =  BDB2 -B2DB  +  7v^(B2D-DB2).f  ev(DB 

=  B(DB  -  BD)B  -  27v777v(DB  -  BD) 

=  B(BD  -  DB)B  -  27v777v(BD  -  DB) 

=  2skw(B2DB-|-27v777vDB), 

where  _ 

£v  =  7v(7v?  —  2IIIy)  =  7v  ^7v7b  +  III^  ,  (7-4) 

and 

B  =  7v"l  -  B .  (7.5) 

^^This  result  follows  from  the  Csyley-Ilamilton  theorem.  An  equivalent  formula,  but  with  more 
complicated  expressions  for  the  coefficients,  had  been  obtained  by  Hoger  and  Carlson  [40]. 

^®Cf.  Hoger  and  Carlson  [40],  Sawyers  [41],  and  Stickforth  [30]. 


-  BD) 

(7.3) 
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The  formulas  (7.1)-(7.5)  also  hold  with  the  replacements  (6.34). 

In  our  derivation  of  (6.33)  we  utilized  the  relation  (6.32)i  for  fi.  By  using  (6.32)2 
instead,  we  find  that  (6.33)  and  (7.1)-(7.5)  hold  under  the  replacements 

(7.6) 

and  that  the  analogs  of  (6.33)  and  (7.1)-(7.5)  obtained  via  the  replacements  (6.34) 
hold  under  the  replacements 

U-U-i,  C-^C-^  (7.7) 

For  a  nonsingular  tensor,  say  A,  let  the  corresponding  letter  in  sans  serif  type 
denote  the  distortional  part  of  A: 

A  ;=  ;  (7.8) 

then  det  A  =  1.  The  distortional  parts^®  F,  U,  V,  B,  and  C  of  F,  U,  V,  B,  and  C 
axe  unaffected  by  the  dilatational  part  of  the  deformation,  that  is,  by  the  value  of 
detF  =  detU  =  det  V.  Note  that  the  polar  decompositions  (1.2),  the  relations  (1.3) 
for  the  Cauchy-Green  tensors,  and  the  formulas  (7.1)-(7.2)  and  their  analog  for  U 
and  C,  also  hold  with  the  replacements 

F_>F,  U^U,  V-^V,  B->B,  C-^C.  (7.9) 

For  any  tensor  H  let 

H„  :=  H  -  1/hI  (7.10) 

denote  the  deviatoric  part  of  H,  so  that  =  0.  Since  the  rate  of  change  of  volume 
per  unit  volume,  or  rate  of  dilatation,  is  given  by 

(detF)7detF  =  4  = /d,  (7.11) 

the  deviatoric  part  Dq  of  the  stretching  tensor  D  is  unaffected  by  the  rate  of  dilatation; 
hence,  Dq  is  a  measure  of  the  rate  of  change  of  shape  or  rate  of  distortion.  Note  that 
the  relations  (1.4),  (6.5),  (6.6),  and  (6.25)  also  hold  with  the  replacements 

D -+ Do,  L  — »  Lq,  La  — >  (l<a)o,  Da-+(Da)o.  (7-12) 

Now  it  is  intuitively  obvious  that  R,  and  thus  D,  should  be  unaffected  by  the  rate  of 
dilatation.  Indeed,  from  (6.33)  we  see  that  the  spherical  part  of  D,  ^41,  cancels  out, 
so  that  D  Ctin  be  replaced  with  Dq.  Likewise,  we  expect  that  fl  should  be  unaffected 
by  the  dilatational  part  of  the  deformation.  Indeed,  on  using  V  =  Illy  V  in  (6.33)  we 

*®The9e  tensors  have  found  useful  applications  in  the  constitutive  theory  of  hyperelastic  ma¬ 
terials  (cf.  Ogden  [42,  Ch.  7),  Rubin  [43],  and  Charrier  et  al.  [44])  and  elastic-plastic  materials 

(cf.  Willis  [45]), 
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see  that  the  Illy  terms  cancel  out,  so  that  V  can  be  replaced  with  V.  In  the  same  way 
we  find  that  the  formulas  (6.26)-(6.33),  (7.1)-(7.5),  (1-11),  (1-12),  and  the  formulas 
for  fl  and  fla,  obtained  via  the  replacements  (6.34),  (7.6),  or  (7.7),  also  hold  with 
the  replacements  (7.9),  or  with  the  replacements  (7.12),  or  with  both  replacements 
together.  In  general,  the  other  formulas  in  Section  6  feiil  to  hold  with  just  one  of  the 
replacements  (7.9)  or  (7.12).  However,  all  of  the  formulas  in  Section  6  as  well  as  the 
formulas  (1.2)-(1.12)  and  (7.1)-(7.5)  hold  when  the  replacements  (7.9)  and  (7.12)  are 
made  together. There  axe  essentially  two  ways  to  show  this.  We  can  either  use 

(7.11)  to  show  that  F  =  LqF,  and  then  proceed  as  in  the  derivation  of  the  original 
results,  that  is,  by  differentiating  the  polar  decompositions  F  =  RU  =  VR.  Or  we 
can  start  from  the  original  results  and  convert  them  by  using  (7.11),  the  relations 

V  = -  I/dV  and  B  = /77b'^'‘B  -  |7dB  ,  (7.13) 

and  their  analogs  for  U  and  C,  which  follow  from  (7.8)  eind  (7.11). 

Since  ■^Do  -  0  and  777u  =  1II\/  =  JIIq  —  IIIq  =  1,  there  is  some  simplification  in 
some  of  the  direct  formulas  which  follow  from  the  replacements  (7.9)  and  (7.12);^® 
if  the  motion  is  isochoric  then  these  simplifications  hold  for  the  original  formulas  as 
well.  We  consider  one  example  here.  Prom  (6.4)  with  the  replacements  (7.9)  and 

(7.12) ,  we  have 

Do  =  Lb[B]  =  -Lb-.[(B-*)‘’].  (7.14) 

These  relations  and  the  direct  formulas  which  follow  from  (5.17)i  were  observed  by 
Leonov  [5].  From  (7.14)i,  the  direct  formula  (5.33)2,  the  identities  (5.40)-(5.43),  eind 
the  fact  that  IIIq  =  0,  we  obtadn  the  following  formulas  for  the  deviatoric  stretching 
tensor  in  terms  of  the  Jaumann  rate  of  the  distortioncd  part  of  B: 

777gDo  =  B(B  -  l^,I)B  +  ,  (7.15) 

where 

=  tr  B  =  tr  B  =  7b  ,  /?2  =  tr(BB)  =  tr  (BB)  =  7b7b  -  77b  .  (7-16) 

The  polar  rate"^^  of  a  tensor  field  A  is  defined  by 

A  :=  A  + An -QA  =  R(RTAR)  RT 

=  A-f  A(n- w)-(n- w)A.  (7.17) 

^^Several  of  these  results  have  been  noted  by  Mehrabadi  and  Nemat.-Na88er  [11]. 

**Cf.  the  formula  for  U  obtained  by  Mehrabadi  and  Nemat-Naaser  [11,  (8.18)]. 

^*Tbcre  arc  a  variety  of  names  u6ed  in  the  mechanics  literature  for  this  invariant  rate.  Here 
we  follow  Dienes  [46],  who  was  motivated  by  the  fact  that  the  rotation  tensor  R  in  the  definition 
n  :=  RR^  arises  from  the  polar  decomposition  of  the  deformation  gradient. 
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In  pairticular, 

V-RURT  and  B  =  RCR'^ .  (V.18) 

o 

By  setting  A  =  V  in  (7.17)3  and  using  the  formulas  for  V  and  f2  —  W  in  Section  6, 

m 

we  can  obtain  a  variety  of  direct  formulas  for  V.  However,  it  is  simpler  to  note  that 
by  (7.18)i,  (6.13),  and  (2.15),  direct  formulas  for  V  in  terms  of  V  and  one  or  more 
of  the  tensors  D,  W,  L,  and  Q  are  given  by  (6.8)-(6.12),  (6.24),  and  (6.38)  with  the 
replacements 


U.-,V,  U-.V,  Dr-.D,  La^L,  W^^W,  (7.19) 

Also,  (7.18)2  and  (6.7)  yield  the  the  simple  formula  B  =  2VDV,  which  was  derived 
by  different  means  by  Dienes  [7]. 

Finally,  we  note  that  direct  formulas  for  the  material  time  derivatives  of  U~^, 
V-i,  V~^,  and  U~^,  and  for  the  Jaumann  and  polar  rates  of  V~*  and  V~i,  follow 
from  the  results  in  this  and  the  previous  section  and  the  identities 

(A-i)- = -A-lAA-^  (A-*)°  = -A-'Aa-*  ,  (A-i)' = -A-*AA-* , 

which  hold  for  any  tensor  field  A. 
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